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We describe fundamental equations which define the topological ground states in the lattice real- 
ization of the SU(2) BF phase. We introduce a new scalar Hamiltonian, based on recent works 
in quantum gravity and topological models, which is different from the plaquette operator. Its 
gauge-theoretical content at the classical level is formulated in terms of spinors. The quantization is 
performed with Schwinger's bosonic operators on the links of the lattice. In the spin network basis, 
the quantum Hamiltonian yields a difference equation based on the spin 1/2. In the simplest case, 
it is identified as a recursion on Wigner 6j-symbols. We also study it in different coherent states 
, representations, and compare with other equations which capture some aspects of this topological 

CNj ■ phase. 

o 

o 

'sj- ■ INTRODUCTION 

A practical guide through the equations of the topological BF phase 

I , Topological order [jy] is a recent key advance in our understanding of the phases of matter. It describes new phases 
H : in which ground states are characterized by the space(time) topology, irrespective to its metric properties, so that 
the physical degrees of freedom are only global. In particular, there is no local order parameter, but new, global 
characteristics, like the fact that the ground state degeneracy is determined by the topology. Topological phases have 
appeared in several condensed matter systems, like those exhibiting fractional quantum Hall effect 0,0] j quantum spin 
^ liquids 0-0] , and also topological insulators pi, loj and topological superconductors [l^, [Hi • Among the new physical 
features, topological order is intimately related to anyons and exotic statistics. 
\ Topological order can be genetically framed into the language of topological gauge field theories describing the 
■ effective behavior, i.e. the low energy regime. Topological field theories are defined without using a metric and do 
not have propagating degrees of freedom. The most well-known is Chern-Simons theory, which is defined in 2 -I- 1 
dimensions. There is another similar theory, known as BF theory, which is defined in arbitrary dimension^, and 
also relies on special properties of the set of flat connections. It is actually quite well understood from the initial 
papers, by Horowitz (l^ and by Blau and Thompson [isj . Further works interested in the symmetries and proving 
the renormalizability can be found in 

[IMS- 
After gauge-fixing, the energy-momentum tensor is in fact BRST-exact. Moreover, the space of non-gauge-equivalent 
solutions is the moduli space of flat connections, whose dimension is a topological invariant. Hence, there does exist a 



X 



?H non-local order parameter determined by Wilson loops which are trivial except around non-contractible cycles in the 
5^ . BF phase. 

The BF model flrst appeared in two dimensions, as the zero coupling limit of Yang-Mills theory [1^. Actually that 
relationship with Yang-Mills theory holds in any dimensions, and it was studied in four dimensions, in the continuum, 
in [l^ (a long time after that phase was observed in lattice gauge theory). The BF TFT appears in relations to 
other interesting field theories. In three dimensions with the appropriate gauge group, it describes three-dimensional 
pure gravity, with or without cosmological constant, making it clear it is exactly soluble |20|]. Even better, the BF 



Plebanski actions [21| offer a way to build four-dimensional general relativity as a BF theory supplemented with 
constraints. The latter are known as simplicity constraints and impose the non-Abelian electric field to come from a 
tetrad. Four-dimensional general relativity again can be formulated a la MacDowell-Mansouri, which is a BF theory 
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^ In 2-1-1 dimensions, it is a Chern-Simons theory. In particular, when the local symmetry group is SO(2, 1), it is a Chern-Simons theory 
based on the Poincare group. 
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for a larger gauge group together with a potential which explicitly breaks down the (topological) gauge symmetries to 
those of gravity j23 |. 

In addition to the effective description of condensed matter systems, and to its relations to the above field theories, 
topological orders provide us with a new way to think of fault-tolerant quantum computations [23'| . The most celebrated 
example is the Kitaev code 24] which is based on a lattice version of the BF model with group IL2 ■ It is also a model 
of importance to test background independent methods in arbitrary dimensions (aiming at a consistent quantization 
of gravity), and in particular the spin network quantization [25l - l27j . 

Spin networks are excitations supported on graphs on the spatial manifold, whose links are colored with spins (or 
representations of a group) and nodes with intertwiners. In the context of Loop quantum gravity, they enable to make 
sense of the quantum geometry of space [HI- They are the same as the string networks of i2y| which describe the 
condensation process leading to topological orders. 

It is known (see for instance [soj). and we will further argue, that spin networks are key objects to get a global view 
on topological aspects which are observed in condensed matter, quantum information and background independent 
methods. For all of them, it is important to have a lattice description, and it is well-known that spin networks are well 
adapted to the BF model since they support an exact lattice version, at least in 2-1-1 dimensions [31U33I]. Furthermore, 
it is a key result in loop quantum gravity that spin network states also span the kinematical Hilbert space of the 
continuum theory when one varies their graphs. Hence, spin networks will be central to our analysis. 

The BF model is often studied with a finite and/or Abelian group (see for a presentation of the background 
independent point of view in the context of finite groups), or a quantum group. Here, we will focus on the (more) 
difficult case of SU(2), which is non- Abelian and in which the Fourier modes are unbounded. In spite of the difficulty, 
we will find that the topological order is encoded and identified by some well-known equations of representation theory, 
which can be handled explicitly by physicists, and which are free from open mathematical issues. 

As the BF theory is a theory of flat connections, the objects which arise naturally in our quantization are re- 
coupling coefficients of the representation of theory of SU(2). These objects, together with spin networks themselves, 
have enjoyed a renewed interest in the last years, in loop quantum gravity, but also in mathematics [35, 36], in quantum 
information [33], and in semi-classical physics (see references in |37l). In particular, one challenge is the understanding 
of their large spin asymptotics. In addition to specific coefficients [381 - 1411] , recent works have unraveled generic methods 
and asymptotic behaviors (i^ - fisj which have nice interpretations in terms of quantized flat simplices. 

In addition to asymptotic approximations, exact results are interesting. A way to capture them is by using recursion 
relations on those re-coupling coefficients (which also provide a way to get the asymptotics [i^, 143] ). As shown in 
[48l - [50| . t hey give a way to probe geometric properties of quantized simplices. In this paper, by extending the seminal 
result of 'si], we show how to get all (well-known, [5l|) recursions on the Wigner 6j-symbol from a Hamiltonian for 
the topological BF phase, and present alternative forms in different bases. At the end of the day, we obtain a practical 
guide through the equations which describe and characterize that topological behavior. 



A new Hamiltonian for the magnetic part 

Our analysis is based on a new lattice Hamiltonian. The usual Hamiltonian [23,[2^[52] consists in two parts, a vertex 
operator and a plaquette operator. The equation for the ground state leads to two conditions: the vertex operator 
imposes the standard gauge invariance, and violations lead to electric charges, the plaquette operator imposes the 
curvature to vanish, and violations correspond to magnetic fluxes. We work with gauge invariant spin networks and 
hence only focus on the second constraint. The new Hamiltonian we introduce imposes flatness of the gauge field 
(hence the topological order) and has the following new features. 

• It is labeled not by a plaquette only but by a pair of a plaquette and a vertex on its boundary. We thus 
have several operators and hence constraints per plaquette. Remarkably, this echoes the proposal of [53] where 
this feature was suggested to avoid the ultra-locality of Thiemann's quantization of the Hamiltonian of general 
relativity in loop quantum gravity. 

• It is built out of bosonic operators acting on nodes, which create, destroy and exchange spins 1/2 between two 
half-lines meeting on a node. 
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• In the spin network representation, the quantum Hamiltonian leads to difference equations on the ground states. 
Here we study them in the case of triangular plaquettes, and find that they are recursion relations known from 
group representation theory. 

• As an operator it can be written in any basis. This gives a variety of equations which all enable to identify the 
topological order in arbitrary dimensions. 

Our proposal is based on recent progress in the context of background independent quantization using spin network 
states. In particular, the spin network formalism has been reformulated using spinor variables 54-60], which have 
unraveled a local U(A^) structure on each A^-valent node. At the quantum level, it is possible to eliminate spins and 
to use spinors instead. We will call the resulting states spinor networks. 

Another progress [3l| concerns the dynamics of the BF model, especially in (2 + 1) dimensions where it also 
describes gravity. The usual way to solve 2+1 gravity emphasizes its topological behavior and makes use of the 
plaquette operator [H, [12, [Ml . However, as a gravitational theory, it can be cast in the form of geometrodynamics, 
with a scalar Hamiltonian which generates time reparametrizations. In that case, it becomes unclear how to solve it. 
We think the plaquette operator is too simple, and solving the model as a gravitational theory is the best way to get 
insights for more complicated models, in particular 3+1 general relativity. 

In (3]| a new Hamiltonian was proposed which mimics as far as possible the scalar Hamiltonian of 2+1 Ricmannian 
general relativity. It is less trivial than the usual projector and plaquette operators since it includes derivative operators 
corresponding to quantum triads. The quantum Hamiltonian equation, i.e. the Wheeler-DeWitt equation can be solved 
in some cases and shown to reproduce the results expected from the topological formulation. The final result is a 
quantization of flat Euclidean geometry. 

The present paper lifts the operator of [31] with the spinor variables to get a new, more fundamental Hamiltonian, 
based on the fundamental representation of SU(2). 

The organization of the paper is as follows. After a short introduction to the BF model in the section HI we present 
in the section [H] the kinematical setting corresponding to SU(2) lattice gauge theory, or loop quantum gravity on 
a graph, in terms of spinors. The dynamics is described in the section IIIII the usual plaquette operator, the new 
Hamiltonian together with its Wheeler-DeWitt equation in the spin network basis, its geometric interpretation and 
the way to recover the usual constraint. In the section ITVl we present the action of the new Hamiltonian on different 
bases of coherent states. We finally compare our key equation with other equations which also encode information on 
the topological phase in the section IVl 

I. INTRODUCTION TO THE BF PHASE 

Let us recall the basics of the Kitaev model [23] for topological order. Every link e of a lattice is attached a spin 
which can be up or down. The Hamiltonian has two parts: star operators (or electric part) acting on nodes s, and 
plaquette operators (magnetic part) Hp acting on plaquettes p, 

eDs eCp 

Hs flips the spins adjacent to a node, while Hp measures the product of spins around a plaquette. One can check that 
[Hs,Hp] = 0. The ground states are those states which are preserved by the all Hg and Hp. On a Riemann surface of 
genus h, the ground state degeneracy is 2^^'') and protected by a gap. 

It is useful to understand the model as a lattice gauge theory for the group Z2. There is an element of Z2 on each 
link. Invariance under Hg means local S2 invariance, which the Gaufi law. Invariance under Hp means that the 
W12 magnetic flux is trivial through each plaquette. This way of thinking generalizes to compact Lie group, such as 
SU(2). Let us put a SU(2) element ge on each link, called its holonomy. Gauge invariance requires invariance under 
SU(2) translations on each node. Triviality of the magnetic fluxes means that the product of holononiies around each 
plaquette is the unit of the group. 



9p=Y[9e = l 

eCp 



(2) 
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The projector onto the solutions writes 



n ^i9p) 



(3) 



plaqucttcs 



where S is the Dirac function over the group (with respect to the Haar measure). It is well-defined for instance on 
Riemann surfaces of genus higher than two and is then independent of the lattice used to define it Il8|: this is the 
meaning of saying it is topological. Conditions for convergence, divergence and invariance are given in [621 ■ 

This lattice model is actually an exact discrete realization of a field theory, called BF topological field theory, and 
known as the zero-coupling limit of Yang-Mills theory [l^ I3| • Though it is a field theory, it only has a finite number 
of degrees of freedom at the classical level which are determined by spacetime topology. At the quantum level, it is 
mostly known from the fact that its partition function is the integral of the analytic torsion over the moduli space of 
flat connections, thus generalizing the two-dimensional case [isj . The continuum action is in dimension d, 



where i? is a su(2)-valued (d— 2)-form, and A a su(2)-valued connection 1-form, whose curvature tensor is F{A). Quite 
clearly, i? is a Lagrange multiplier imposing the vanishing of F{A). In addition to the standard SU(2) gauge invariance, 
there is another gauge invariance (which actually contains diffeomorphism invariance), in which B transforms like 
B 1-^ B + dAfj. This symmetry is responsible for the disappearance of all local degrees of freedom. The path inte gral 
treatment was done in fisj , and most of the Hamiltonian analysis and canonical quantization in [l2j . We refer to [8l4ll| 
for the relevance of this theory in the effective description of topological insulators and topological superconductors. 

Thanks to its simplicity, it provides spin network states with a dynamics we are able to control [12, [E^i in the case 
d = 3 0. The situation is more subtle in higher dimensions and under current investigations. And in spite of its 
apparent simplicity, it produces an interesting connection with a non-trivial equation from Lie group representation, 
the Biedenharn-Elliott identity. 

The relation to the SU(2) lattice model is the following. The holonomies ge along the links are the Wilson lines of 
the gauge field A. The SU(2) action at each node of the lattice comes from the effect of gauge transformations on those 
Wilson lines: they transform it only on its end points. Hence, the Gaufi law indeed asks for translation invariance at 
each node. The field strength F{A) is regularized as usual in lattice gauge theory, by Wilson loops. In particular, the 
flatness equation F{A) — becomes the statement that Wilson loops around the plaquettes are trivial. 



The classical, kinematical setting we use is just the phase space of lattice gauge theory. It actually coincides with 
that inherited from loop quantum gravity restricted to a single graph. The formulation we present was proposed in 
(S?! . Further classical and quantum kinematical aspects can be found in [s^ - IsgI [58l - l60| . 

Let r be a closed graph with L oriented links and V nodes. We introduce a classical phase space formed by attaching 
a spinor € (D^ to the source vertex of each link, and another spinor ie G (C^ to the target vertex of each link. The 
Poisson brackets read: 



often short-handed to ±. To make sure we construct SU(2) invariant quantities, we define the following, standard 




(4) 



II. KINEMATICS OF LATTICE SU(2) GAUGE THEORY 



A. The lattice spinor variables 




^ In three dimensions, the situation is well controlled for open manifolds, or topologies of the type a Riemann surface cross a real interval. 
But it is still unclear for generic closed manifolds because the gauge-fixing of !32] then leaves a residual, non-compact gauge symmetry. 
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contractions of spinors. First define some notation: 

k) = (ll) , (^1 = an* i^+r) , (6) 



together witli the duality map 



\z]^\,z)=^^l_^i j, i.e. (,z)- = (-f)^--(z--)*. (7) 

The natural invariant contraction is; (wiz), and using the duality map we get a second invariant contraction: 

{w\z) =Y,iw^)* z^, and [w\z) ^J2{-l)^'^w-^ z^. (8) 

A A 

One checks easily that: {gw\gz) = {w\z), for any g G SU(2). The invariance of [w\z) is ensured by the fact that: 
{(;gz)'^ — gAsi'iz)^ , which means that <;z transforms like z under 811(2)0. 

Forming scalars with these inner products enables to define invariant observables, but we are also interested in 
covariant objects. Let us consider vectors. We need obviously to tensor two spinors, and then use an intertwiner 
Hi ~^ to map covariantly the tensor product to the representation of spin 1. Such an intertwiner is 

well-known and its components are basically the matrix elements of the Pauli matrices. More precisely, take the 
matrix- valued vector a which reads in the fundamental representation: ax = ( i o ) > ''"y ~ ( i ) ' ''^^ ~ ( o -i ) • The 
commutation relations are: [cr*,cr-'] — 2ie*"'j,(T''. Then, vectors, known as the flux variables, are given by: 



{z\a^\z)^-\ (9) 

A,B 

and they satisfy the 5u(2) algebra. 

More than vectors, we can actually reconstruct the classical phase space T* SU(2)^ from the spinors (and this is 
the lattice gauge theory phase space) . The idea is to introduce a group element g^ to each link which sends the spinor 
on the source to that on the target, 

W) \Ze) ^ , \Ze] \Ze] ,-,„^ 

ffe = ~ _ :, and ffe = = — F===. (10) 

\/{Ze\Ze) ^J{Ze\Ze) \J ^Z^\z^) {Ze\Ze) 

The spinors Zg , actually completely determine ge [STj , 

\Ze) {Ze\ + \Ze] [Ze\ . 

9e = _ _ . (11) 

^J {Ze\Ze) {Ze\Ze) 

Those objects are the holonomies (they can be seen as parallel transport operators, from the source to the target, 
arising as Wilson lines of a gauge field). The full underlying symplectic geometry has been described in [57] (though 
with different conventions) . The Poisson brackets are: 

{XlXi)=^\Xl {Xlg,)^g,T\ (12) 

where the matrices (Ti)i=i_2,3 are anti-Hermitian generators satisfying [t*,t^] = -^t^ . All other brackets vanish. 
Hence, fluxes X\ act as left invariant derivatives, and should be thought as attached to the half-link given by e and 
its source vertex. 



^ The reader may have noticed that this second contraction uses the antisymmetric e-tensor in the fundamental representation, e = ( J ) . 

Its matrix elements in the irreducible representation of spin j are: (i, Tt) = (— l)^~"<5_„,m , and they will be used in a systematic 

way in the remaining of the paper to contract Wigner 3jm-symbols. 
^ They read: r* = —icr*, in terms of the Pauli matrices (c'). 
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Right invariant derivatives are obtained by transporting to the target vertex of e via the adjoint action: 

Xe = Ad(ge) Xe - .9e Xe ff"', (13) 

together with: 

{Xl,Xi} - -e'\X^, {XI, g,} = g,. (14) 

Notice that the adjoint representation on the algebra is exactly the representation of spin 1, acting on 3- vectors. 

That holonomy-flux algebra is the starting point of lattice gauge theory. In this context, one introduces a SU(2)^- 
action, by group translation on each node v of F. The Gaui3 law generates those transformations through the Poisson 
brackets, and local SU(2) invariance is imposed by the constraint 

X,- J2 ^^ = ^- (15) 

outgoing eDv ingoing eDv 

Straightforward quantization on T* SU(2)^ leads to the Hilbert space i^(SU(2)^), equipped with the Haar measure 
(or L^(SU(2)^/ SU(2)^) when gauge invariance is required). This is the kinematical Hilbert space of loop quantum 
gravity restricted to the graph T in 3+1 dimensions (see [l^, [l^,!!^ for reviews, or more specifically [1^), but also in 
2+1 dimensions f32l. [52|. 

The phase space T* SU(2)^ assigns six real variables to each link, while the spinor setting gives eight variables. 
The relation between both frameworks is given in j57| . The idea is to perform symplectic reduction with respect to 
a U(l)^-group action. As spinors and their quantization are well known, it is very appealing to reformulate (at least 
the kinematics of) loop quantum gravity this way, as done in [s^, and hopefully get new insights on the dynamics. 

A key difference between spinors and holonomies/fiuxes is that the latter encode the information on the links of 
the graph, while spinors enable to factorize it on the nodes. Therefore, we now focus on a single A^-valent node v 
and consider without loss of generality that all the links which meet there are outgoing. Following [s^ , we introduce 
elementary SU(2)-invariant local observables (i.e. acting on the node). We use here latin indices a,b,c,d to denote 
the legs of the node. 

Eab = {Za\zb), Eba = E^b, (16) 

Fab = {Za\zb], Eba = -Fab (17) 

That set of observables forms a closed algebra. Interestingly, the observables Eab satisfy a u(iV) sub-algebra, and this 
is the reason why this formalism has been coined the U(7V) formalism for loop quantum gravity [5J-|56|, 

{Eab,Ecd} — — i{5bcEad — SadEcb), 

{Eab,Ecd} ^ - i{Sbc Fad - SbdEac), 

{Eab, Ecd} ~ 0, 

{Eab,Ecd} — — i{SacEdb + Sbd Eca — Sbc Eda — SadEcb)- (18) 



B. Quantization: Intertwiners and Schwinger's boson operators 

Let Hj denote the carrier space of the irreducible representation of SU(2) with spin j € N/2, of dimension dj = 2j+l. 
From the Peter- Weil theorem we known that L'^{S\J{2)) = (Bj^^Uj ® H*. 

The Hilbert space on F is L'^{S\J{2)^), and thus it is spanned by the products of the matrix elements of the 
holonomies in any set of L irreducible representations (je) attached to the links of F. Therefore, a matrix element of 
the holonomy ge along e in a representation ke acts as an operator on the Hilbert space by a simple multiplication, 

{ke,p\ge\ke,q){je,n\ge\je,m) = {ke,p\ge\ke,q) {je,n\ge\je,m) 
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V 




FIG. 1. The tetrahedral graph we consider throughout the paper. The orientations are the same as those of the graphical 
representation of the 6j-symbol [sj. The three fat hnes form the cycle (126) we will consider in order to explain the action of 
the new Hamiltonian. 



The second line just comes from Clebsch-Gordan re-coupling, to map the reducible tensor product ® to the 
sum of irreducible representations (Bj''^y^ The quantities into brackets are Wigner 3jm-symbols (equivalent 

to Clebsch-Gordan coefficients but more convenient for our purposes) (sij . 

The flux variables become insertions of anti-hermitian generators (left or right derivatives) , 

K {je,n\ge\je,ni) = {je,n\geT'- \ je,m), XI {je,n\ge\je,'m) = {je,n\T' ge\je,m), (20) 

with — — ^cr*. 

A basis of the gauge invariant Hilbert space on F, _L^(SU(2)^/ SU(2)^) is given by spin network functions. The 
latter are formed by contracting the magnetic indices of the Wigner matrices along each link with a specific tensor at 
each node. This tensor has to ensure the invariance of the function under (left or right) translation of the holonomies 
and it must then be an intertwiner between all representations meeting on the node. Consider typically a 3-valent 
node, where the links are outgoing and carry the spins ji, j2jj3- There is a single intertwiner Lj-^j^j^, i.e. a single 
invariant vector in Tij-^ (g) Tij^ (8) Hj^, which is up to normalization the Wigner 3jm-symbol[3. It is convenient to see it 
as a map and write its tensor elements like: 

{ji,rm;j2,m2;h,rm\L,,,,,jO} = (^^^ , (22) 

where |0) is the normalized vector in the trivial representation. 

If some links are ingoing, as it happens on the node where (62,66,64) meet in the figure [U one dualizes the 
representations, here J2, J4, and writeqj 

04,m4|.,j,.,,b2,m2;j6,m6) = (-l):'-™^ (-1)^-™. (^J^^ J^^ . (24) 

Gauge invariant observables are obtained from traces of holonomies, and more generally from full, rotation invariant 
contractions of flux indices and magnetic indices of Wigner matrices (see [3l| for a recent use of such operators). 

Since the flux carries a vector index, any gauge invariant observables containing it will add a link to F in the 
representation of spin 1. This operation is known as a grasping and enables to extract quantum geometric information 



^ Being an invariant vector in Hj^ ® "Hjj ® means 

E (il,-lMii,-i>02,n2|,b2,™2>(i3,n3|9b3,rn3>(;^^^ t) = (ni Z il) ' ^^'^ 

mi ,m-2 I'fi'Z 

for any g e SU(2). 

^ Note that with these conventions the standard group averaging formula for the tensor product of three representations is: 

I dg 55?=! D'^'-Hg) = 'jij2J3lO)(Okn.j2.M* • (23) 

where D^j^{g) is the Wigner matrix. 
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from spin networks [3l[. Taking tensor product, one easily creates graspings with arbitrary integral spins. However, 
we would like to perform graspings with half-integers also. This is what lifting the phase space to spinors allows. Let 
us proceed to their quantization. 

The spinor we have introduced on each half-leg of F is obviously the classical version of the Schwinger representation 
of the angular momentum using two harmonic oscillators. Upon quantization, the latter become two annihilation 
operators {o.'^)a=-i/2,+i/2: such that: 

[a^,a^t] ^^AB .^^ [a^,a^]=0. (25) 
From them, the su(2) generators can be constructed as quadratic operators: 

r^ = -^{a\a^\a). (26) 

These objects are anti-hermitian: (t*)^ = — t*, they indeed generate the expected algebra [t^, Tj] — e^^^ Tk, for k — 
x,y,z. (a'^) can be thought as an operator- valued spinor which transforms under the fundamental representation of 
SU(2). 

A natural basis (|?i-^)) is that obtained by diagonalizing the occupation number operators, (a^)ta^ and (a+)^a+. 
This basis is exactly the standard basis of spins and magnetic numbers (|j, m)) on the sum over all irreducible 
representations Qj^^Hj. The relation between both sets of eigenvalues is: 

j==i(n+ + n-), m=i(n+-n-). (27) 

Hence, we can write the action of {a^,a^^) on the standard basis. Compared to the generators (r^) of the algebra, 
the new feature is that they enable to go up and down between different irreducible representations. Indeed, 

I J, m) = ^J] + (-1)3 --4^ \j_}_^m-A), a^^ \j,m) = j + {-l)^-^m + l\j + ]^,m + A) . (28) 

It is well-known that the matrix elements of the generators (r^) can be re-expressed as Wigner 3jm-symbols with 
a spin 1 (that is the spin carried by the generators). This is the way one could construct a Hamiltonian operator 
generating a recursion relation on the 6j-symbol, well-known in re-coupling theory, in (sH ]. Actually, a similar result 
holds in the present case: 

{k,n\ a^Uj,m) =v/d-d^(-l)^-"+3-A+i S,^^^^ ^ _Q , (29) 

with the notation dj = 2j ' + 1. It also suggests that their action can be translated into re-coupling equations. 

Now we equip the graph F with these operators. At the quantum level, to take into account SU(2) invariance (the 
Gaufi law), the degrees of freedom at each iV-valent node are intertwiners. However, as the Schwinger's bosons change 
the spin, we cannot consider the spin on each leg to be fixed. Instead we take the invariant tensors living on N copies 
of ©jeN^j, 

7^^ = 0Inv[^f^iH,J , (30) 

{Je} 

where one sums over all sets of irreducible representations attached to legs. We have assumed that all legs are oriented 
outward^. Schwinger boson operators enable to define operators on that space which were out of reach in the usual 



Changing the orientation of a link amounts to changing the corresponding spin j with its dual representation, denoted j* 
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loop quantum gravity framework. The operators we are interested in are quadratic and create, annihilate or exchange 
quanta between two given legs of a node: 

A=±k 

Fee' = [ae | ^e' ) = J] (-l)3-^a;^ ® fl^ , F,, , = -F^e' ■ (31) 

These operators are clearly the quantization of the observables Fee' and Fee' , (fT6|) , (jl7p . Due to the special contractions 
of the spinor index, they are invariant under the global SU(2) transformations acting on the node, i.e. they commute 
with the generator J = X]e=i^e- The set of operators F, F, F^ form a closed algebra and like in the classical case, 
the operators Fee' satisfy a u(iV) algebra. 

Putting the Schwinger operators on spin networks will just mean for us that we have that structure on each node, 
acting on Hn- One link e has two sets of operators, {af) on the source vertex, and (af) on the target vertex of e. 

As we see in the definition ([5]), the symplectic structure on the incoming legs has the opposite sign to that of 
outcoming legs. This is reflected in the quantization by the fact that and a^^ act on the right, i.e. on bras. This 
is due to the dualization of the representation. Thus, our conventions imply that those operators satisfy the usual 
algebra: 

[a^,a^^] =S^^ id . (32) 

Hence, we can forget the tilde when considering their matrix elements. But the interpretation is reversed: since they 
act on the right, the creation operator, i.e. that shifting the spin on the leg by +i, is a'^, while a,-^^ is the annihilator. 



C. Graspings with spin 1/2 and holonomy operator 

We show here that the operators F, F, F^ generalize the usual operators used in the context of loop quantum gravity, 
by producing graspings with the spin 1/2, and we collect the results of their action. 

One way to evaluate these operators is by using the standard action of the creators and annihilators A brief 

examination shows that one has to use some recursion relations on the 3jm-symbols to carry out the calculation. 
Another way through, which instead highlights the re-coupling nature of those operators and their interpretation as 
graspings, can be used thanks to (|29l) . From the results below, closed forms are obtained by taking the values of the 
re-coupling coefficients in a textbook [51], which are reported in appendix for completeness. 

Let us write the action of, say, the operator F21 on the node where ei, 62, 63 meet in the figure [1] 

E21 Cnj.jM = 4\ji,mi)^ af\j2,m2) ^ \j3,ms) (^^ ^] , (33) 

74,7711 ,7712 

The last line is obtained after some re-coupling. The operator F12 thus performs a grasping on the intertwiner tj^j^j^ 
with a spin 1/2 inserted between the spins ji and j2- The shift of ±^ on ji, J2 comes from selecting the value fci in 



First rewrite the action like 
-^21 tjiisjslO) 

where the object l'=i'=2.'3 is 

7-'=lfc2J3 _ 



ni ,712 



A,rni ,7712 



-A nil 



ki 
-ni 



i-iy 



k2 j2 \ ( jl j2 j'i 



(35) 



(36) 



We have used first the definition, and in the second equality the expression of the Schwinger's operators in terms of 3jm-symbols I I29II . 

represf 

( nj n2 ) ■ '^oi'^ precisely, the pre-factor comes from: 
(fci , m ; A:2 , n2 ; is , "13 1 tfei fc2 j3 |0) • (37) 



2^nin2m3 is a tensor with three indices, ni,n2,m3, respectively in the representations of spin A:i,fc2,i3. It is easy to check that it is 
SU(2) invariant, and hence proportional to the 3jm-symbol, In}n2m3 



ki k2 is 
Pi P2 P3 



E7-fcl'=233 
P1P2P3 

P1,P2.P3 

That pre-factor is built with four 3jm-symbols, which are contracted following the pattern of the tetrahedron, and thus produces a 
6j-symbol. That 6j-symbol is special since it has a fixed spin being 1/2. The exact evaluation yields the desired result. 
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FIG. 2. The 3-valent nodes represent 3jm-symbols, where legs carry the spins. A magnetic index is summed when there is a 
link joining two nodes. The action of E21 is a grasping between ei and 62, denoted by a dashed line which carries the spin 1/2. 
The final result is proportional to the 3jm-symbol with the spins fci = ji — i, A:2 = 72 + |. 



the tensor product: 

1^,,®^! =^,,_i + (38) 

and the value k2 in Hj^ 18 Hi. The values of ki, k2 are determined by the operators, here ai and a^- The contraction 
of the magnetic indices A, mi, 7712 results in an invariant tensor in Hki ^ 'Hk2 Hj^, which is obviously proportional 
to the normalized intertwiner ikik2ja- The pre- factor involves the quantity into curly brackets which appears in the 
final form and known as a Wigner 6j-symbol with a spin fixed to 1/2. This is summarized in the figure [21 

This calculation produces a spin 1/2 grasping which generalizes the standard spin 1 grasping, well-known in loop 
quantum gravity. The spin 1 grasping simply comes from the action of the fluxes. Instead of considering E21 — {z2\zi), 
one can just consider the dot product of the corresponding fluxes, Xi ■ X2. The quantization of this observable is 
obtained by promoting each X' to an anti-hermitian generator acting on the node. The calculation has been 
reported precisely in [3l|, and leads to 

(n • r2) i,,,2,,\0) = (-1)^'^+^'=+^'^+Vji(ji + l)d,, + l)d,, 1^:^ ^,,,.,3|0). (39) 

The fact that no spin gets shifted in this case is clear from the fact that the operator (ti • T2) selects in the tensor 
product "Hji iSiHi the factor Jij^ . The operator (ti • T2) is quadratic in the generators and thus quartic in the operators 
a"^. Hence, it can also be obtained from some square of the spinor contraction £'12- A similar process was detailed 
in [3l| to get graspings with a spin 2 from those with the spin 1. It obviously involves the square of the spin 1 
operator together with some symmetrization. At the quantum level, it can be derived using the recursion formula 
on the 6j-symbol to relate the different sets of eigenvalues and pre-factors, i.e. here to express { ^-^ } in terms of 

k2 j-2 ja ' ' 

Let us compute in addition the action of a operator on the same node, — (ai|a2], 

F2I hmjsl^) = «f ® (-I)^"^«^^^|j2,m2) ® |J3,™3) 

A, mi ,7712 ,m3 

^ (-l)n+,^+^^+i^d,,dk,d,2dk2 + [ll 1} ^k.k2,M- (40) 

The spins ji, J2 are both raised to ki — ji + ^,k2 = j2 + ^ due to the creators. 

Quite clearly, all operators E, F, F'^ will generate a 6j-symbol with a spin 1 /2, of the form | / | > with ke — je^\- 
The main difficulty is however to establish properly the overall signs, which are so important to get the correct 
recursion formula in the end. Let us here report the results of a few actions, in particular on nodes with incoming 



Jl J2 J3 
mi 7712 fTT-S 
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and/or outcoming legs. The conventions and notations correspond to the situation of the figure [TJ 



'2 J6 is 



'-fclfe2j3|0) : 


(41) 


'-fclfc2j3|0) 7 


(42) 




(43) 




(44) 




(45) 




(46) 




(47) 




(48) 



EiehUih = (ail«6) Lj-^f^je = {-^y^^'^'Wdjidk, djjk^ 6k,,n-k^k<,^o<,-^ j^^ J^j 

It is essential that all usual loop quantum gravity operators on a fixed graph can be expressed using the operators 
E, F, F^. We have already discussed the case of flux operators acting on a single node, as quartic operators in terms 
of spinors. Since the loop quantum gravity phase space {Xe,ge} can be recast into a (constrained) phase space based 
on spinors, one should also be able to quantize the holonomies (ge) with the Schwinger's boson operators. 

On the kinematical Hilbert space Z/2(SU(2)^) over F, a matrix element of any holonomy in a any irreducible 
representation acts by simple multiplication. We consider the elementary situation: the holonomy in the fundamental 
representation along a single link e carrying the spin je- The multiplication is re-expanded onto the matrix elements 
of irreducible representations as follows, 

(l/2,A|5e|l/2,B)0e,n|ge|je,m) 

= E dj^i-^y--"^ (1 (-1)'='"' (I „m) (Je,iV|5e|Je,A0 ■ (49) 

N,M 

There are only two terms in the sum, Jg = je ± i, and for each of them the 3jm-symbols can be understood as matrix 
elements of Schwinger boson operators thanks to (1^ . 

{l/2,A\g,\l/2,B){j,,n\g,\j,,m) = ^ \{j„n\af g,af^ \je,m) + {j„n\ (-1)5" 5e (-l)^~''ae"'' lie,m)l . (50) 

This may be written in a more compact, symbolic form, 

geD^-^^Hge) = (|Se)(ae| + |Se][ae|) -^D^^^^ge). (51) 

This is obviously the quantization of the classical expression of the holonomy in terms of spinors (llip . Notice that the 
factor 1/dj which enters here is nothing but the quantum translation of the holonomy normalization 1/-^/ {'z\'z){z\z). 
However, one has ordering ambiguities when promoting the classical expression (|lip to an operator. What we find 
here is that the result of the correct ordering is 1 /dj . 
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The relation ([?T|) enables to get the translation of the parallel transport of spinors (ITUl) at the quantum level: 

^{9e)AB {je ~ l/2,n\ge a^\je,m) = (je - 1/2,77,10-^ ge\je,m), 
B 

Y.^9e)AB {jeMde (" 1) ^ fl-^t | _ 1/2, „,) = (^^ , 1) I ^-^t |^-^ _ 1/2, ™) . (52) 

B 

The analogous formula for the spin 1 is equivalent to the well-known fact that the adjoint action on a Pauli matrix 
is a rotation on the matrix- valued vector f. Indeed from: g^^r'^g = J2j -^iyVj'''^ ' '^^ S^^- ^(5)*j 0' "-Is ~ 
{j,n\T' g\j,m). 



III. DYNAMICS OF BF THEORY 

The previous section presents the kinematics of loop quantum gravity on a single graph, or equivalently the frame- 
work of lattice gauge theory, formulated with spinors. We now study the dynamics of the BF model. The Hamiltonian 
consists in two constraints, one imposing the Gaufi law, which we have treated as part of the kinematics, and the other 
the vanishing of the curvature, i.e. of magnetic fluxes. 

First, we present the direct equivalent to the plaquette operator of the Kitaev model, which is based on characters 
of the Wilson loops. Those calculations are not really new, but the geometric interpretation as tent moves is and was 
suggested to us by B. Dittrich. Then we summarize the ideas coming from the new Hamitlonian introduced in [Slj . 
before presenting its extension to spinors. 



A. The plaquette operator: Tent move evolution 

We consider a plaquette p bounded with n links. The constraint gp ~ 1, has three real components, since 
dim SU(2) = 3, on each plaquette. This can also be seen by gauge fixing all holonomies around p to the unit but one, 
which is then set to 1 by the constraint. 

The set of constraints is sufficient and the projector (1551) is well-defined for example on Riemann surfaces of 
genus higher than two [18]. To implement it on the lattice, it is standard to expand it over the SU(2) modes of the 
delta function, 

^(ffp) = '^ip Xjpigp) , (53) 

where Xj is the character (the trace) in the spin j, normalized to Xj(l) — dj. All terms in the above sum are well- 
defined as SU(2) invariant operators on spin networks, and this is the way Noui and Perez solved 2-1-1 gravity in 

However, the set of constraints on the Wilson loops usually has redundancies which make the projector ill-defined, 
in the absence of a suitable regularization, as it would involve divergent products of Dirac deltas. One can still 
look for physical states directly as solutions of a well-defined Wheeler-DeWitt equation. As suggested by the above 
decomposition, one can use the character of the holonomy gp around p, for any irreducible representation of spin j, to 
enforce the constraint 

N 

Xjidp) IV')phys - dj |V')phys = , *^ y ■ (54) 

When written in the spin network basis, this constraint leads to a difference equation (sol . [6^ which has some 
interesting features. The matrix elements of each holonomy in the character can be re-coupled with those of the same 
holonomy in the spin network function. This shifts the spins of the cycle (ji, . . . , j„) to (fci, . . . , fc„) with sums over 
ki = \ji — j\, . . . ,ji + J. One is left with graspings in the spin j at each node of the cycle, and they can be factorized 
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FIG. 3. A pictorial representation of (|55p . The character Xj along the closed loop acts on the left. On the right we have 
depicted the situation after re-coupling. A specific 6j-symbol is extracted on each node, and one has to sum over the colorings 
k\, . . . ,kn. The dashed lines correspond to the dual 2d triangulation to the plaquette if we think of the latter as embedded in 
flat 3-space. The vertex s of the 2d triangulation is then dual to the plaquette. 




FIG. 4. Here we have displayed the geometric interpretation of the character operator on the plaquette as a tent move. The 
vertex s is evolved to a new vertex s', with an edge of length j + the tent pole. Between the initial and the final triangulations 
we have a piece of 3d triangulation. The character operator then generates the evaluation of the Ponzano-Regge amplitude on 
this triangulation. 



from the function as 6j-synibols with a spin j. Precisely, this gives 

ki,...,kn e=l ^ ■' ^ ■' ^ ■' 

This is represented in the figure [3] Note that this enables to vifrite the projector on the topological ground state 
restricted to a single plaquette in the spin network basis, 

j,fci,...,fc„ e=l ^ ^ ^ ^ ^ ^ 

(56) 
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It has a very nice geometric interpretation in terms of a (2 + l)-dimensional evolution. The reason for the (2 + 1) 
dimensions is that it is natural and always possible to consider a plaquette to be dual to a piece of a 2-dimensional 
triangulation embedded in flat 3-space. The plaquette itself is dual to a vertex s, and the links on its boundary are 
dual to edges which connect s to vertices (si) dual to the surrounding plaquettes. 



Then we use the Ponzano-Regge interpretation [SJl of the 6j-symbol: it is a weight associated a tetrahedron whose 
edge lengths are given by (je + 1/2). A Ponzano-Regge amplitude is then obtained by taking products of 6j-symbols 
corresponding to a gluing of tetrahedra, and summing over the spins of the internal dual edges, keeping those on the 
boundary fixed. 

This interpretation shows that the equation (1551) is actually computing the Ponzano-Regge amplitude for a piece of 
a 3d triangulation built as follows. Draw from s an edge outside of the 2d triangulation, with length {j -1-1/2) going to 
a new vertex s* , and connect s* to the vertices (s.^). This process creates a piece of 3d triangulation with n tetrahedra, 
all sharing the edge (ss*). Such an evolution process of the canonical surface is known as a tent move, and it is 
depicted in the figure |4l It gives an evolution in discrete time steps, where the triangulations of the constant time 
slices are all the same. This evolution has been analyzed in background independent approaches to lattice gravity in 
[65j |. The tent pole is here (ss*) and has length (j -I- 1/2). By Fourier transforming ([55l) . one sees that the evolved 
state Tp{ke) is obtained by summing over all admissible values of the spins on the initial surface, 

il,-jTi e=l ^ ^ ^ ^ ^ ^ 

In particular, a physical solution to the flatness constraint must satisfy 

(Tent, ^A) (fee) =dj V(fce) , (58) 

that is it must be invariant under tent moves weighted by their Ponzano-Regge amplitudes. A way to enforce this 
requirement is to use the projector on the topological sector (|56|) which is geometrically interpreted as a summation 
over all possible tent pole lengths. 

B. The spin 1 Hamiltonian 

However, there are some drawbacks. In particular, the difference equation one gets shifts the spins on each link in 
the boundary of p. This means that one cannot solve the Wheeler-DeWitt equation on each single face independently 
of each other. And yet we know that in some situations, one should be able to get and solve an equation on each spin 
independently. 

The simplest example is that of the canonical surface S'^ triangulated by the boundary of a tetrahedron. Let us 
write explicitly the spin network function, labelled by six spins (ji, . . . , jg), following the figure [TJ 



mi , . . .,7716 

ni,...,nQ 





h 


i, / 




TO2 





„{jc}f„^ _ I 31 32 33 \ I Jl J5 36 \ I 33 34: J5 W 32 J6 34: 

-ni -ris me) I-tt-s -n^ m^J 1-^2 -rig 



n( 1 



{3e,ne\ge\3e,rne 



(59) 



The range of summation is —j^ < me,ne < je for each link. The natural inner product on i^(SU(2)^/ SU(2)^) is 
evaluated with the Haar measure on SU(2)^, and given the chosen normalization for spin network functions we get 
here 



' 5 



(4c^t^i4c\=^>=n^- (60) 

Given the trivial topology on the 2-sphere and the chosen fixed graph, there is a single 'physical', i.e. topological, 
state, which is peaked on fiat connections with full measure: 

tpphysigi,---,96) = S{96g594)S{g6gig2^)S{g^^g4g2)- (61) 
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Here 5{g) is the Dirac delta on SU(2) with respect to the Haar measure. To get its expansion on the spin network 
basis, we take the inner product with a spin network function sl^^^{gi, . . . , ge), 

V'phys(jl, . . . ,J6) = / n '^^'^ 4ot^(ffl> • • • >56) '0phys(5l, • • -,96), (62) 

e=l 

■ (63) 

J4 J5 J6 J 

This resuh simply comes from the fact that the physical state enforces the evaluation of the integral on the flat 
connections (up to gauge transformations), while there: s^j^y{\) = {6je}- 

Thus we could call Wheeler-DeWitt equation in that situation any equation which characterizes the 6j-symbol. As 
a re-coupling coefficient, the latter is known to satisfy a special identity, the Biedenharn-EUiott, or pentagon, identity 
(5l| . From it, one derives (46j a recursion formula which enables the evaluation of the 6j-symbol, 

The coefficients are given in (sij for example, and they are built on special 6j-symbols which have one spin equal to 
1. It turns out that the coefficient A+i (respectively A-\) vanishes when the spin ji reaches the highest admissible 
value (respectively the lowest admissible value). This means that the recursion is initially first order, i.e. it can be 
iterated from a single initial condition. 

The relationship to the topological, flatness constraint ([2]) was unraveled in [3l|. We want an equation acting an 
a single face which selects the 6j-symbol in the case of a cycle with three links. It was shown in [3l| that the above 
recursion can be derived from the quantization of a constraint on the phase space of loop quantum gravity. To build 
the constraint the idea is to extract the matrix elements of the holonomy around a plaquette g^ by projecting onto a 
flux variable on the left and on the right. Consider the situation depicted in the figure [TJ with g^ = g^gxg^^ based on 
the vertex where 62 and eg meet. We obtain a spin 1 Hamiltonian by projecting the rotation matrix R{gp) onto the 
fluxes X2,Xq, 

Hille, = E ^6 " ^(56 5i 52"')^j) M = X, ■ X2 - ■ kd{ge 31 g^^)X2. (65) 

Ad is the adjoint representation on fluxes seen as Lie algebra variables, which is of course equivalent to the standard 
vector representation of SU(2). Thus, we look at the rotation R{gp) in the basis spanned by the flux variables. The 
constraint can be re-written in a more convenient way using the parallel transport relation between left and right 
fluxes (1131): 

<'le. =X^-X2~X^- Ad(5i) X2. (66) 

On a spin network function, the first term Xq ■ X2 is diagonal and given by (l39l) . The second term is a bit more 
involved but an essential point is that there is only one holonomy gi, and hence only one spin, ji gets shifted by 
— 1,0,-1-1. This leads to the difference equation (sij : 

^+i(ji) ^(ji + 1) + ^o(ji) V'(ji) + ^-i(ji) ^(ji - 1) = 0. (67) 

The coefficients are: 

AoUi) = (-1)J"^+J'*+J« -^2 ). + (_i)2ji+j2+j3+j5+j.(2ji + 1) <j I \ { ■'I I \ , (68) 



^±i(j-i) = (-i)^^^+^"^+^^+^"^+^"«+^(2(j-i±i)+i) ^'f^^ ]^ I 



ji±l ji 1 [ jji±l 31 1 
je 36 35 



(69) 



This difference equation is exactly the recursion relation satisfied by the 6j-symbol. 
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FIG. 5. A graphical representation of the action of the new Hamihonian. The basic idea is that the holonomy around a closed 
loop in the topological sector only depends on its homotopy type, so that we can deform the grasping on the left to that on the 
right, picking up this way some holonomy which must be trivial. 



A graphical representation of the action of H^^'' is given in the figure [H The idea is that on the topological sector, 
the holonomy along a closed loop only depends on its homotopy type. Hence, we can make a grasping, i.e. insert a 
line inside the plaquette whose ends are glued on the boundary of the plaquette, and being in the topological sector 
means that the curve can be deformed, since the plaquette is homotopically trivial. The Hamiltonian proceeds like 
this: it inserts a infinitesimal curve (i.e. with trivial holonomy) at node between two links (62,66 on the figure [5]). 
The constraint is just that this curve can be deformed to run all along the plaquette, picking up the corresponding 
holonomies. Hence, it indeed produces a constraint on these holonomies, whose product has to be trivial. 

In addition to this topological interpretation, it enjoys the following properties. 

• It is labeled by a cycle (which bounds p) and a vertex v in the cycle. Since a cycle has at least three links, this gives 
a sufhcient number of constraints to impose the topological condition R{gp) = 1 G S0(3) (i.e. gp — ±1 e SU(2)). 

appears as a regularization of the scalar Hamiltonian in the 2+1 gravity, 

Hon = EfE^^\F^, . (70) 

F is the field strength, and E the pull-back of the triad to the canonical surface, and a, b are space indices. 
It comes from projecting the constraint F — onto the lapse and the shift. This splits the full 3d diffeo- 
morphism algebra into a constraint which generates spatial diffeomorphisms, plus i?GR which generates time 
reparametrization. However, it is unclear how the same splitting can be done on the lattice. 

The answer provided by H^^^ is to project the curvature around a node, say s in the figure [31 onto the normals 
to each of the triangles sharing s. Generically, if there are at least three triangles, all directions of 3-space will 
be spanned. This way we get an evolution of the vertex s in at least three "time" directions. 

• In 3d, the phase space of loop quantum gravity consists in Regge geometries, which means they have intrinsic 
geometries determined by the spins interpreted as lengths. The constraint Hi]p — imposes a relation between 
the intrinsic, spatial geometry and its extrinsic curvature characterizing the embedding into flat spacetime. 

• In 4d BF theory, the geometric interpretation is similar in the Regge sector, but they are also non-Regge, 
discontinuous intrinsic geometries in the phase space |66| . known as twisted geometries. Hi]^ provides a well- 
defined constraint and well-defined recursion relations for the twisted geometries [1^. 

• The asymptotics of the equation can be studied, for example, using coherent intertwiners. The solutions to 
the lowest order in the WKB approximation in the Regge sector are linear combinations of exp(±i5'R), where 
S'r is the Regge action of the simplex. We thus get a criterion to tell whether the Regge action is part of 
the asymptotics of a lattice model [i^. The Hamiltonian Hi]l is an example which generates exactly solvable 
recursions beyond this asymptotic behavior. 

Those are features we would like to reproduce using spinors, for the group SU(2) instead of S0(3). 
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C. The spin 1/2 Hamiltonian 

1. The classical spmor Hamiltonian 

We labeled the operator h'^'I in with a superscript (1) because the derivatives transform in the adjoint 
representation of SU(2), that is with a spin 1. As a consequence (see [3l[), the recursion relation obtained involves 
shifts on a spin by the amount ±1. 

Thus, if the recursion (j67|) is iterated from an initial condition with all spins being integers, we can solve the family 
of 6j-symbols with integer spins. This amounts to solving the model for S0(3). However, from this solution we cannot 
get information on the family of symbols with, say, three half-integer spins and the three others being integers. Hence, 
these two families of 6j-symbols have to be related by hand, choosing appropriately (i.e. "fine-tuning") the initial 
conditions. 

If one could derive a recursion with spins 1/2, one would then be in position of solving directly the model for the 
group SU(2). To this aim, one needs variables transforming under the fundamental representation of SU(2), in contrast 
with the flux variables. We have just described such variables in the previous section, the spinors. 

Hence, we mimic the construction of the constraint (|65p . In particular, the graphical representation should be 
unchanged, as given in the figure [SJ but with insertions of spin 1/2 instead of spin 1. Technically, the rotation matrix 
R{g) should now be changed with the fundamental representation. Thus, we propose 

H^eee'X = C^dl^^] 1^6 1 96 9192^ " 1 \z2) ■ (71) 

Note that we could actually write four such constraint, by changing 'z2,'zq with <;'z2 and/or ^zg, 

1^2], 

\Z2). (72) 

1^2]. 

Some simple algebra on spinors gives that the sum (with correct signs) of these four observables leads back to the 
holonomy around the cycle (126) of our tetrahedral graph, 

Hi-;e\l + Hi+-l - Hi++l - Hi-^-l = (?2|2-2)(Z6|Z6) tri(.g6gi52-i - 1) , (73) 

where tr^ is the trace on 2x2 matrices. Therefore, our observables provide a decomposition of the holonomy 

operator around the closed loop (126) into simpler components, which will translate into quantum operators generating 
specific shifts of the spins around the loop. 

Moreover, these observables are not independent and there are some relations between them. For instance, taking 

into account that = — (z|g|z'] for group elements g G SU(2), one gets ijigliel = He^eie^'^^\ for 62,66 = ±- 

Furthermore the fact that g = gegig2^ is unitary, gg^ = 1, translates into quadratic relations between the iJ*^^^) 
observablesQ- Finally, we obtain that there are only three independent real components. More precisely, we would like 
to impose all the constraints = 0. These are 2x4 real constraints. Due to the complex conjugation relations, 

this already reduces to 4 real constraints. Then discarding the pre-factors, (26122] and so on, imposing = 

amounts to requiring that the four matrix elements of the 2x2 matrix gi g^^ ~ 1 vanish. Since ge gi is in SU(2) 
and only depends on three independent real parameters, we do have in the end only 3 independent real constraints, 
which simply express that the holonomy around the closed loop (126) must be the identity, gegig2^ = 1- 

The above form of the constraint contains the holonomies along three different links. But it can be simplified as 
depending of the holonomy along a single link with the following lemma. 



We decompose the relation gg^ = 1 for g = 96 9i 92 ^ explicit matrix elements by inserting twice the identity in term of spinors 

{z\z) 1 = 1^) (2I + 1^] [^1 for z = Z(; and z = Z2. 



-f^ee'etes = (^e\^V [zel 96 9l 92 " 1 
-f^i^eVi = [^6 122] (2^61 56 5132^^-1 
-ffi^eteo = [^6|22) (26| 56 5l52"^ " 1 
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Lemma 1. The observable Hegeie2 admits a form with only the holonomy along ei, 



Hi;-l = {ze\z2] J ff'lf^l^^l [ze\ 9i ^2) - {ze\z2] [ze\z2) , (74) 

V \Z2\Z2) [ZqIZq) 



and a form completely in terms of spinors, 

I {z2\l2){ze\l,) ^,,, J^6|(|^l)(^l| + |?l][^l|)k2) 

They simply come from applying the parallel transport relations of spinors, (fTOj) and (fTTj). to get rid of the different 
holonomies. 



2. The quantum spinor Hamiltonian 



Promoting the above quantities to operators is done by using the quantization map we have discussed in the previous 
section. Mainly, the spinor contractions are sent to operators Eee' , Fee' and F^^, . 

At the quantum level, one faces numerous ordering ambiguities. However, several of them are resolved from the fact 
that we know the standard action of matrix elements of holonomies on spin networks. From this, we have shown on 
the equation (jSlI) that the norm of the spinors in the expression of 51 , 1/ -y/ (zil^i) (zi|zi), only contributes to a factor 
1/djj in the spin network basis. In the equation (j52p . we have shown that the parallel transport [zq\ does not bring 
any factor, which means that the ratio \/ {z2\z2) /{z2\z2) is naturally set to 1 at the quantum level. 

This way, one has reduced the ordering ambiguities to a single one, the ordering of the two factors in the definition 
((7T|) of the Hamiltonian. Let us define two operators, one for each ordering, say left and right. Using the forms of the 
classical observable displayed in the lemma [1] the quantization with the left ordering is 

^l^eV 4it^ = ([ael 31 |«2) - PelS^)) 4ct ^ , (76) 

(Fei E,2 + E,e F12) - ^62) sl^^ , (77) 



while the right ordering is 



^LillR ^t^ot ^ = (["el 51 l«2) - [ae\a2)) {ae\a2] s^t^ , (78) 
= (Fei ^12 + ^16 i?i2) - -F62) F^, sli;^ . (79) 

Note that i^|g commutes with (Fgi £'12 + -B16 -F12), since i^2g acts on the vertex where 62, eg meet, while that (Fgi E12 + 
Eie F12) act on the two other vertices of the cycle. Hence, the true ambiguity is the ordering between F^g and Fq2- 

The quantization of for a, /3 = ± is obtained by changing in the above definitions 0,2 (and/or ag) to <;a2 

(and/or c^ag). Then, it turns out there is a relationship between the two orderings, 

f^("'«,J^=i7(-"--^^ (80) 

V 666162111 / 666ie2|L \ ' 



This is the quantum equivalent to He^eie2 — -ffe6e"e2 

Theorem 1. Consider the 2-sphere triangulated by the boundary of a tetrahedron, where kinematical states take the 
form ipigi, ... ,55) = [IlLi '^jj V'(ji, ■ • ■ , je) Stet^(ffi' ■ • -^de)- The Wheeler-DeWitt equations for the physical, 

topological state are 

Hi:eZlLW=0^ (81) 
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for a, P ~ ±1/2. In the spin representation, they give difference equations, 



^(Ji + 1/2, J2 - /3, . . . , J6 - a) + - 1/2, J2 - /3, . . . , J6 - a) 

+ V'(ji,-.-,j6)=0, (82) 



which are known to come from the Biedenharn-EUiott identity and determine the 6j-symbol. 



The coefficients of the difference equations are 



J2 J2 + ^ J3 



2 i 

je J6 + a is 



(83) 



The relation (j80| between iJ|L and the adjoint of H^Yi, translates to identities between those coefficients, 

Ai"'^\ji,j2 -a,.. ~(3) = (-1)"+'^ ^-"'-^^(ji, j2, • • ■ ,J6) , 

^,,+1 A^^fUi + ^, J2 - a,je - /?) = d,, 4T'"^^(ji, J2, . . . , Je) , (84) 

which can also be checked directly from ([83]) . Those identities can be used to prove ([82| by direct calculations, which 
are reported in the appendix IbI 

There are four possibilities for (a,/3), and hence four difference equations for each plaquette. However, due to the 
discrete symmetries of the tetrahedral graph, only three equations are actually different. The physical reason why 
three equations per plaquette are independent is that the Hamiltonian is built to impose (classically) the equation 
gp = 1. For each face, this is a SU(2)-valued equation which has three real components. 

Combining three of these equations, one can reproduce "4^ the spin 1 recursion formula (j64p on the 6j-symbol. Since 
the latter can be solved from a single initial condition, we can obtain families of 6j-symbols where the spins differ from 
the initial spins by integers. Two such families are \ "^^ ^ \ and i "^^^ "2+2 I .^^j^gj-g ni,n2,nQ G N. 

To relate the different families, one has to resolve steps of half-integers (which are not integers). This is precisely done 
using the spin 1/2 recursions, which consistently set the relationships between the initial conditions of the different 
families. 

We will prove the theorem for H^^^\ as it works similarly for the other cases. It follows from the two propositions 
below, together with the equation ([5(I| . 

Proposition 1. ^^g~e7e2|R '^tct^ annihilates the tetrahedral spin network on the identity, 

(V'phy. k^>>> = sli^Hg,, . . . ,56)|,.=i = . (85) 

Before showing up the calculations, let us emphasize that it is easy to understand the action of the operator. 
Indeed, consider the form (|78l) of the operator. The multiplication of the spin network by matrix elements of gi in its 
fundamental representation induces shifts of the spin ji to ji ±1/2. Then the role of the operators ag, ae, ^2, 02 is i) to 
make the whole operator SU(2) invariant, and ii) to shift j2,j6 so that the rules of tensor products of representations 
are still satisfied after the shifts of ji (i.e. it is necessary to shift j2 so that J3 is in the tensor product 'Hj^j^i®'Hj^_i)- 

The proof of this proposition is based on a very simple method which can be applied in more generic situations than 
the tetrahedral spin network [67| . 

Proof of proposition [TJ First observe that the effect of F^^ — (aQ\a2], given by the equation (|48p . is to shift the 
spins j2, ie on the half-lines of 62, cq meeting at v to I2 — J2 — \ — J6 ~ \- Thus, the intertwiner ij^j'j^ in the spin 
network function is changed with iffj^. Up to some multiplicative pre- factors, we are left with the remaining action 
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of the operator, that is [ag|gi [02) — [a6|a2), usmg the form of the operator given in (175)) . It is not hard to see that 
this gives zero due to the somewhat trivial following identity, 



J2 i-l)^-^{l/2,A\gi\l/2,B) {I2,n2\g2a''\32,m2) {k.ne\gea~^ \je,me) 



A,B=±i, 



(86) 



We emphasize that this is one of the key identities in this paper. It states that we can commute the creators and 
annihilators with the group elements, here 52,36, when the latter are taken to be the unit of the group. Moreover, in 
that case, the sum of the left hand side simplifies since (1/2, A\ gi |l/2, = 6a, b- 

To see the action of [agl^i |a2) — [061^2) in the above equation, one has to dress it with the matrix elements 
(je, ^e|<?e|je, 'Tie)]! On the links 6 ^ 62,65, and contract them with the four intertwiners tjifajs,'- 



H*i*j4- That gives the equivalent form 



J- X! [^-^6 ^ l/2,"-6|^61 hiJ'Jeljl + l/2,n-i;j5,"-5) {jl + 1/2,^1; ^2 - 1/2, ^2; J3, "-3|£^12 tjuajs |0) 
•^^ ni,...,n6 

+ Ue - l/2,?T.6|£'i6 tjji-jelji - l/2,ni; j5,n5) {ji - l/2,ni; j2 - 1/2, ri2; ja, '^3|-Fi2 tjusjalO) 
X (J5,^^5|^J*J•J5b3,?^3; ^4,^4) (j4, n4|t(j,_ 1 ).jjj2 - 1/2,^2; - ^,ne) 

[Oe, W6kj-jjj8|ji,mi; j5,m5) {ji,rni;j2,m2; j3,m3\Lj^j^j^\0) (js, m5|ij.j-j, jja, ma; J4, m4) 

nil T.-ime 

X (i4,TO4|^62i(j2-i)*06-i)*i4b2,™2;j6,™6) • (87) 

Both sides of that identity are contractions of four intertwiners following the pattern of the tetrahedral graph. The 
holonomics have been evaluated on the unit. The left hand side contains two different terms, obtained by putting 
for the holonomy gi = {\ae){ae\ + \ae][ae\) / dj-^ in the above (|86)l. These terms correspond to [agl^i \a2) — (-Ei6^i2 + 
F6iEi2)/dj-^, and encode two shifts of ji, by ±1/2. The right hand side is the part Fg2 = [a6|a2). □ 

Proposition 2. Moreover, expanding the different terms of the proposition [T] using re-coupling produces a spin 1/2 
recursion formula on the 6j-symbol, 



J4 



J5 J6~ 2 



Jl J2 J3 
k js J6 



+ a'^T\ji) 



jl -\ j2-\ ja 



J4 



]5 J6 



Proof of proposition [2j It simply follows from considering the form of the operator h'^^J^^^^^-^ displayed in the 
equation ([79ll , and evaluating it using the actions of the operators E^e' , F^e' , F^^, which are given in the section III CI 
Using (^7)) . we get that Fq2 -PJe diagonal on spin networks. 



^62 J^26 *tct 



j2 j2-\ \ 

je - h je ji 



(89) 



The term which lowers ji by -1/2 is obtained thanks to (gB]), (gl]) and 

±E,eF,2F^, 4>> = i-iy^+'^'-''i-iY^+'^'-'^i-iy^+'^+'^d^,^.d,,d^^_.^ 



^31 
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The term which increases ji by +1/2 is obtained using P5|) . (|42p and (jlS]) again, 

^F,,E,,F^, = (-l)^-^+^-^+^«+i(-ir^+-'>-+^'H-lF^+^'^+^'^d,, + i rf^^^ 



-^31 



(91) 



Then, we sum the three contributions, add a phase (— l)J2+j4+i6+i^ divide by a common factor, dj^dj^_i djgdj^_i x 
{ I ''^ . ^ ^ \, and finally evaluate the spin network function on the unit of the group. □ 



D. Geometric interpretation: Quantum Euclidean geometry 



Although we consider throughout this paper the group SU(2) as the local symmetry group, our methods are expected 
to generically apply to any compact Lie group or finite group (see [34| for a recent review on spin net and spin foam 
models with finite groups). But SU(2) is interesting from the geometric point of view. In the three-dimensional context, 
it is the (universal cover of the) isometry group of tangent spaces to spacetime. This leads to an interpretation of the 
spin network geometry which is well-known in the context of 3d quantum gravity, revisited here with the new spin 1 
and spin 1/2 Hamiltonians. 

The key idea is to consider the piece of 2d triangulation dual to the spin network graph, whose triangles are dual to 
the 3-valent vertices, edges dual to links, and vertices dual plaquettes. By a trivial extension of the graph (with links 
carrying the spin zero), we can assume that it is a triangulation of the canonical surface (space), hence characterizing 
its topology. 

The usual phase space variables, holonomies and fluxes, characterize the geometry of an embedding of the trian- 
gulation into 3-space. The Hamiltonian constraint of curvature vanishing asks for an embedding in fiat 3-space. We 
briefiy sum up how this can be seen and refer to [3l| for details and references. 

The three fluxes on a vertex are interpreted as the normals to the edges of the dual triangulation. This is possible 
due to the Gaufi law ([T5l) . which is exactly the closure condition for the triangle dual to the node. That also means 
that the norm of a flux, X^, is the squared length of the dual edge. Similarly, the dihedral angle (j)ee' between two 
dual edges (e*, e *) of a triangle is contained in the dot product of Xe, Xg', where the dual edges e, e' meet at a vertex: 

COS^ee' = -Cee' ' '^^^^ 

where Cee' is 1 if e and e' are both outgoing or ingoing, and —1 else. Lengths (and dihedral angles), which are 
determined by fluxes, completely characterize some intrinsic geometry of the canonical surface, and are invariant 
under local rotations changing the frames in which fluxes are expressed. 

From the 2d angles (and hence lengths), one can compute the dihedral angles between neighboring triangles as if 
they were embedded in flat 3-space. Consider triangles ti,i2 dual to vertices wi,W2 of the spin network graph, which 
are linked by the edge e. The classical angle of flat space between the two triangles ti,i2 can be computed from the 
three 2d angles around the node which is dual to the cycle c, using the formula: 

O (V V\ COS^eiea - COS^eie COS^ese .„„s 
cos etit2 (A , X) = — — . (93) 

Sm0eie Sm0e2e 

Further, the holonomies ge also captures some notion of dihedral angles, through the formula 

a ^ ^ • Ad(.ge)A^t2 

cose.,,AX,X,g)^-^^^-^j^, (94) 

where Nt are normals to the triangles which are evaluated like N^^ — e* ^.X^X^^ . 
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It turns out the Hamiltonian Hi]l simply expresses the difference between the two notions of dihedral angles between 
neighboring triangles, 

=sin(^eie sin^e^e (cosetit2 - COS^t^f^) ■ (95) 

Hence, the topological equation for the ground state, Hi]l = 0, just imposes that the dihedral angles 9t-^t2 computed 
from the holonomies are those of the flat embedding. 

Quantizing if^^^ hence means quantizing flat Euclidean geometry. This is what the Wheeler-DeWitt equation (p7|) 
does. 

At the quantum level, the re-coupling coefficients yl_i.o,i thus provides a well-defined way to evaluate the quantum 
operators which correspond to the above products of trigonometric functions. (This is particularly clear in the large 
spin limit [H, S^). The most transparent case is that of the 2d dihedral angle, e.g. between the dual edges to 62,65 
in the figure [T] We have seen that the spin network function is an eigenvector of X2 ■ , so that 

/ {ie}N/ N + 1) +i6(j6 + 1) - jiik + 1)] {jj , , 

(cos^e^ee ^tot ) (ffi , ■ • • , ffe) = ^ , . , . . , . 4 (9i, ■ ■ ■ , 9e) ■ (96) 

2VJ2(j2 + 1) JeUe + 1) 

This exactly reproduces the classical expression of the dihedral angles of a triangle, with quantized lengths £e{je) — 

It is reasonable to think that the spin 1/2 Hamiltonian carries a similar geometric interpretation, but with half 
dihedral angles instead. It is however hard to make precise, since one needs operators diagonal on spin network 
functions, so at least quadratic quantities like i^62^26' ^^'^ then one faces ordering ambiguities. So the best we can do 
is to check that the large spin limit is indeed the same as that of coscjjeiee/'^ computed from the above operator. 
The latter gives 

^^g2 4'e2ee COS 062 66 + 1 



^ {Vhij2 + 1) + VmUg + 1) + Vuiu + 1)) WJ2U2 + 1) + VjeUe + 1) - VuUi + 1)) 

4v/j2(j2 + l)j6(j6 + l) 

on spin network functions. From the spin 1/2 graspings, we get 

F62 _ {j2 + J6 + ]i + 1) (j2 + k - k) 



^22^66 (2i2)(2j6) 



(98) 



In the large spin limit, the square roots of Casimirs go to the spins themselves, j{j + 1) ^ j, and the two quantities 
above have the same limit. 



E. Equivalence with the flatness constraint 



We would like to see exactly how the recursion can be used to find wave functions satisfying the fiatness everywhere. 
As it was already observed in [s^], it is clear that the recursion relations we have written so far only depend on the 
fact that we consider a cycle with three links. In particular, they do not depend on the fact that we have closed the 
graph by gluing the links 63, 64, 65 to a node in the figure [T] 

Hence, we now consider a cycle with three links and 3-valent vertices, and use the notations of the figure [TJ but 
assuming that the links 63,64,65 do not necessarily meet. We still have exactly the same recursion relations, one 
for each node of the cycle. This gives three independent difference equations to impose three real constraints on the 
holonomy around the cycle. With their help, we want to show that the physical state factorizes on the cycle. In 
the spin network basis, the state is a function of the spins. In this configuration, its dependence on the spins ji, J2, Je 
actually completely factorizes: 

V'phys(je) = (-1)^-'= g} 0(j:), (99) 
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where the spins jg are the spins of all links of the graph except ei, e2, eg. This is known in the 3d gravity case from the 
projector onto physical states [s^ . Such a factorization property is actually very natural. Indeed, a flat connection 
induces on the boundary of a face with three links holonomies which are up to gauge necessarily trivial. Then, it is 
known that spin networks evaluated on the identity satisfy such factorizations, and the case of cycles with more links 
is also known (50| . 

The idea to extract the factorization is that the recursion relations (for spin 1/2 and hence spin 1) hold on ji, j2, Je. 
We can thus implement them from V'(jij J2, J6, Je) so as to reach an initial state with ji = 0, j2 = js and je = js. The 
factorization ([99)) is then proved with the initial conditiorF^ 4'ije)\ji=a = (^1)^""^ | jl | 0(je)- 

Now we are ready to present the consequences of the factorization (|99p in the group representation of wave functions, 
in which we also get a factorization of the cycle, 

'0phys(ffl,g2,ff3,g4,55,56,...) = S{g6 gi §2^) 0(.93ff^\54,5655,...) , (100) 

with 

^9392 ^,94,9695,- ■ ■) = X![n^J^] '^(•^e) S^''^^ (53^2" \ 34, 9655, ■•• ) ■ (101) 

We thus see that the goal has been achieved: the wave function has support on holonomies with trivial SU(2) parallel 
transport on the cycle, 969i92^ — 1- 

To obtain the above equations, one writes the group Fourier transform of the physical topological state, 

V'phys(ffe) = [ll'^Jo] ^Phys(je) S^^^HSe) , 

{>} 

= T.[l[d4 Hfe) E 1^:^ (-l)^-.^»(ffe). (102) 

{je} e' jlj2j6 

Using the group averaging identity J dh D^^^'^ (h) (g) D^^^^h) (g) D^^^^ (h) = i-ja2h\^)i^\''hhj3, and the Fourier expansion 
of the Dirac delta on the group, S{g) — J2j '^jXjig), the sums over ji,j2,j6 can be explicitly performed, 

E 'ind.Ae {^-^ f\ (-1)'''^ s{>>(5e) = 6{ge9i 92^) ■s^''^H9392 \ 9i, 9695, ■■■)■ (103) 

ji -.32 -Je 



IV. NEW REPRESENTATIONS OF THE TOPOLOGICAL EQUATION 

We are now interested in finding alternative forms of our key equation ([82)) . In the usual spin network basis, it is a 
difference equation which corresponds to the well-known recursion formulae for the 6j-symbol. Beyond the result of 
having a new Haniiltonian, with its nice geometric interpretation, one of the major progress is that the key equation 
comes from an operator. Hence one can simply re-write H\4>) = in a different basis, and the physical content, i.e. 
the flatness constraint of the topological model, is preserved. 

One of the technical interesting aspects of using the Schwinger's bosonic representation of SU(2) is that it is possible 
to go to bases of coherent states easily. We will first recast the recursion relation in a basis of SU(2) intertwiner, closely 
related to those introduced in |68] . and which are detailed in [s^ ]. This basis is mixed, i.e. it contains spins and spinors. 
Then, we will move to a basis based on the coherent states of the Schwinger's bosons, thus getting rid of spins. The 
constraint equation will appear as a holomorphic differential equation solved by Schwinger's generating function of 
6j-symbols. 

Further representations are provided in the appendix [C) 



It will be convenient in tlie following to explicitly factor a phase (—1)^^^. The basic reason is that in our conventions es is oriented 
inwards the cycle while 62,64 are outwards. 
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A. The recursion formula in the coherent SU(2) basis 

For z = ^ ^+ ^ e (D^, introduce the following states, 

\J,^) = ^ 10) = 7(27)1 E ,1 , b-'-) • (104) 

They are normalized to {j,z\j,z) = (jz-P + jz+P)^^. They are just the standard coherent states for each of the 
two quantum harmonic oscillators, projected onto the subspace of fixed spin j. In particular, they are holomorphic, 
homogeneous polynomials of degree (2j). The annihilators and creators act as 

a^\j,z)^^z^\j~^,z), and a^t |j- ^ _i_ ^ 1^- + 1 , ^) . (iqs) 

In addition to lowering the spin by 1/2, the annihilator multiplies the state by z^, while the creator a"^^ increases 
the spin by 1/2 and derive the state with respect to z^. 

The duality map ([7]), allows to define anti- holomorphic states, \j,z] = |j, <;z). We will mainly use their duals on 
which the boson operators act as 

[j,z\a^=^-^=^ ^[j + lz\, and [j, z\a^^ ^ ^ {^1)^-^ z'^ [j - ^, z\ . (106) 

We construct a basis of intertwiners by group averaging tensor products of these states, naturally coined coherent 
intertwiners, 

H>)(^e)|0) = / (g) 5|je,^e). (107) 

isU(2) V 

and from them, we define the spin network functions in this basis. On the tetrahedral graph of the figure [I] they are 
labeled by six spins (je) like before, but also one spinor z^ for each source node s(e) of each link e, and another spinor 
Ze for the target node t{e). Consider for example the node of the figure [1] where 61,62,63 meet, all three outgoing. 
Then the intcrtwiner components are proportional to the 3jm-symbol, 

{Ji:m^;n,m,;js,ms\cu ^)\0)=('' T IT v^^^^C^iS^^ f'' , (108) 



1,2,3 ' 



and the proportionality coefficient is an invariant holomorphic polynomial on (21,22,23). For additional interesting 
aspects on these intertwiners (especially of higher valence) , we refer to [i^ [t^] ■ The polynomial is actually exactly 
known [5l|, and is the generating function for the 3jm-symbols with fixed spins. That gives 



'-iU2i3 (^11^2,2:3) — -fjlj2j3 (^1:^2,^3) '■j2j2j3 (109) 



with 



3 



P,,,,,3 (21,22,^3) = ... . . ,]j n yM:y- Nkl)^-^-'^'^ [^3k2)^-^-'^'^ [Z1|Z3)'^-^-'^'^ , 

^Ulj2j3j (.Jl23 + Ij! 

— (110) 

^^j^j^^^) = / (-^123 - 2il)! (Jl23 - 2j2)! ( Jl23 " 2j3)! 



(Jl23 + 1)! 

We have set J123 = ji + j2 + js for the total spin. On a vertex with some incoming links, where target spinors live, 
it is convenient to built the coherent intcrtwiner using the bras [je, Ze|- For example, on the node where 61, 65, 65 meet 
on the figure [TJ we take 

''jrj5i6(^i'^5,Z6) = / dh [ji,zi\h~'^ <^[j5,Z5\h~^ <^h\je,ze) . (Ill) 
Jsv{2) 
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This convention gives formula similar to the one above, 

Equipped with these intertwiners, we can form the corresponding spin network function on the tetrahedral graph, 

4 6 



ide) = rr TTbeitWel /lt/e^5e^s(e) |je,^;e) , (113) 

= [P]i]2]Azi,Z2,z:i)Pj^j^^j^{zi,z^,ZQ)Pj^j^j^XI:^ s\ll^ {gi, . . . , g^i) . (114) 

This factorization between the standard spin network function, which is independent on the spinor labels, and the 
holoniorphic polynomials, which are independent of the arguments of the function but capture the whole dependence 
on the spinors, is specific to the case of three-valent nodes. In particular, the evaluation on the trivial holonomies, 
5e = 1, gives the 6j-symbol times the four polynomials, which we denote like 



J4 % i^^^^'d = [Pnj233 izi,Z2,Z3)Pj,j,j^{zi,Z5,ze)Pj.,j^j,{z^ j-^^^ -^^^ -^^^j . (115) 



We are now ready to state the main result of the section. 
Proposition 3. The equation ^^e6eie2|-R'^tor^^'^^"'^(ffe)|i = reads 



2ji + 1 



'-A,B=±l/2 



+ (2j2)(2j6)(2ji) [zi|z6)[zi|z2) ^\r5 ^2-^ .j;3 



- (2ji + 1) ^ ^) f f {z^rze) = . (116) 

The remarkable feature of this equation is that it obviously generates the same shifts on the 6j-symbol as in the 
previous recursion, but the coefficients Ao^±i/2 of the recursion have been exchanged with simple multiplications and 
derivatives with respect to the spinors. 

Proof of the proposition [Sj First, it is clear from the argument already used in the proof of the proposition [T] 
that h'i „l \ T}s\iV'^'"'^''\9e)\i is indeed identically zero. The key relation is displayed in (j86| . 



666162 l-R tct 

To actually evaluate the action of the Hamiltonian on the coherent spin network, we use the form (17^ , which we 
recall: (^(Fei E12 + Ei^ F12) — ^52) F^^. It is thus sufficient to compute the actions of the elementary operators 

E, F, F'' in the basis of coherent intertwiners. Notice that these operators are invariant under SU(2), so that they 
commute with the group averaging in the definition of the coherent intertwiner (|107p . For example, 

^26 hSJSM (^2, 26, Z4) = / dhV] (-1)^-^ [j2, Z2\a''^'' ® [jg, ^eja^t /j-i ^ _ (II7) 

^SU(2) ^ 

Then, we just apply the creators and annihilators in this basis, to get 

^26 HiilJi (^2, Z(i,Zi) = V2j2 2j6 \z2\ze) (i2, Sg, Z4) ■ (118) 

Remember that i^Je ^1^^ ?a:si operator in the right ordering of the Hamiltonian. We can actually factor (and 
ignore) the multiplicative contribution, ■\/2j2"2jg [z2|26), and consider that is only here to shift the spins J2, J6 to 
j2 — 1/2, j6 ^ 1/2. Indeed, by dressing the trivial identity (|86)) with intertwiners, we derived the equation ((87|) . Both 
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and ([57]) hold in any basis, and in particular we can take the coherent intertwiners in ([57]) . Hence, we can just 
apply the latter, and plug into it the following actions of the other relevant operators, 



Fl2 Z2, Z:i) = yj2ji 2j2 [Zi\z2) ij,_ij^_ij^izi,Z2,Z3) 

-^12 ''J1J2J3(^1' ^2, Z3) 



2j2 Vf— 

1 + 1 ^\dz^ 



d 



2j 



A Z2 j ij 



i + ^i2-ij3(^l'^2,2:3) , 



^61 Lj'^j;je(zi,Z5,ze) 



i{zi,Z5,Z6) 



E16 h'ljfjeizi, Z5, ze) = [zi\z6) t(j-^_i).^-.jg_i(zi,Z5,Z6) 



(119) 
(120) 

(121) 

(122) 
(123) 

□ 

Another proof using the special factorization (|109p of the three-valent coherent intertwiners is available, though it 
is unclear how it could be more generally applied. The idea is that the coefficients of the recursion are taken into 
account as multiplications and derivatives with respect to spinors. But the dependence of the coherent 6j-symbol on 
the spinors is only through the holomorphic polynomials P. Hence, the latter satisfy the key relations which enables 
to get the correct coefficients. Let us just give an example. We have seen just above that Fq2 produces a second order 



Fe2 i-u^^^r Ue -^^rjA^2,ze,Zi) 



dze 



dzo 



(^2,^6,^4) • 



operator, [■§^\-§^)- Then, it is tedious but straightforward to show that 



1 



V2jj2j^ 



d 



dze 



d 

dz2 



, (z2, Z6, za) = y (■?2 + J6 + 3a + 1) (j2 + je - ji) Pj^^y^^y, (z2, zq, Zi) , (124) 
= -4"'"^ ^dj,dj^^idjjj^_i Pj^_ij^_ij^{z2,Z6,Z4) . (125) 



The coefficient Aq is given in ([55)) . All other terms can be evaluated this way, and the pre-factors, in particular 
Pj2-^jg-^j4^ can be factorized, as well as the other polynomials, finally leading to the spin 1/2 recursion relation on 
the 6j-symbol. 



B. The Schwinger's generating function for 6j-symbols 

In this section, we will get rid of the spins, and fully transform the recursion relation to a differential equation with 
respect to the spinors only. On each node of the spin network graph, we build some intertwiners in the following way. 
Assume that the node has its three links outgoing, like the links 1, 2, 3 on the figure [1] Consider 

Liz,,Z2,z,) ^ ti' / ^ ^("'^l^^) 10)'^ ' (126) 

ilj2j3 lle=l(2je)- J e=l 

- (Jl +j2+j3 + 1)! . ^ 

— 2^ 1 3 ■ ~ — hij2j3[^1^^2,Z3) , [LZ/) 



n:=i(2je 



where the intertwiners Lj-^j^j^{zi, Z2, z^) are defined in (|107p . We thus get invariant vectors on (8)e=i {®jJ~LjS)^ the 
tensor product of the spaces of all irreducible representations on each leg. When some lines are actually ingoing on 
the node, we use the same convention as in the previous section, e.g. on the node of the figure [T] where 2, 6, 4 meet 

f~ ~ . ^ (j2 +i4 +J6 + 1)! f~ ~ ^ .,„o^ 

t(z2,Z6,^4)= 2^ . . tj.j|j,(z2,Z6,24) ■ (128) 
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On a spin network graph, we associate one such intertwiner to each node, and glue them using the holonomies on 
the hnks. This is the natural gluing, where all irreducible representations are orthogonal. This leads to the notion of 
spinor networks, which we define as 



^ nLi(2je)! 



i9e) 



(129) 



where the coherent spin network s^^ 
node V, 



Ze,ze} jg defined in (|113p and Jy is the sum of the three spins meeting on the 

Jl23 = jl + j2 + js , Jl56 = jl + jb + je , >/264 = j2 + je + J4 , -^345 = js + ^4 + js • (130) 

Its evaluation on the unit is a 12-spinor symbol, which is the Bargmann representation of the 6j-symbol, mostly 
known as the Schwinger's generating function of 6j-symbols, 



n 



(J. + l)! 





J2 




IJ4 


J5 


ie J 



nni 



\Jv-2j, 



(131) 



The reason why it generates the standard 6j-symbols is the special factorization of the coherent three- valent intertwiners 
(I109p . There are twelve brackets [ze"\ze') in the final product over vertices and links. The notation is as follows. To 
each vertex v we consider a cyclic ordering of the links which meet there, respectively (123), (156), (264) and (345j3- 
In the above formula, we have written (e, e', e") the links meeting at each vertex with the correct ordering. Since there 
are three choices of e, we have three brackets for each node, and since there are four nodes, that gives a total product 
of twelve brackets (and the tilde is understood for the spinor on the target vertex of each link) . 

This generating function is known to have a closed form [5l| . which can be derived from several methods, a la 



Bargmann (see [7l| for example) using complex integrals, or a la Schwinger 72[ using his boson operators to evaluate 
inner products of states, 

'^^^^'^"^"n AM -M^ ' -^(^e,^e)= J2 I[i^e'\Ze)+ Ul^-'l^e), (132) 

[l+A4[Ze,Ze)\ 3-cyclcs v 4-cyclcs v 

with 

E [z2\zi)[zi\ze)[ze\z2) + [Z3\z2)[z2\z4)[z4\z3) + [zi|z3)[z3|z5)[z5|zi) + [z4\ze)[zQ\z5)[z5\z4) , 

3 — cycles v 

E ^ [^lk3)[24|23)[z4|z6)[zi|z6) + [Z3\Z2)[Z3\Z5)[Z6\Z5)[Z6\Z2) + [^^2 l^l) [^2 1^4) [^5 |24) [zs l^l) ■ 

4 — cycles V 

There are two types of contributions, one from the cycles of three links, and the other from cycles with four links. 
There are obviously four cycles with three links, corresponding to the triangles on the tetrahedral graph. There are also 
three cycles with four links. Each of them is obtained by removing a couple of two opposite edges on the tetrahedral 
graph. That gives the cycles (61626465), (61636466), (62636565) on the figure[TJ On each node v of each such cycle, the 
function SlzejZe) picks up a contribution from the bracket [z'^\Zf,) where e,6' are the two links meeting at v on the 
boundary of the cycle. 

Proposition 4. The equation H^^^^^^j^S^^ 
generating function, 

d_ 

^dzf 



'(5e 



is true and corresponds to a differential equation on the 



[ E (l?-^)( 



A,B=±l/2 



9zf 



9. E E^-^flzc) f/9^J/9zn) (^+E^5;^5^a) 



eDi)264 C 



[z.\z,)[z4\z2) (2+i ^ E-?|c)(2 + ^ E E 



dz^' J idze\dz2 



eDV26i C 



eD'Ui23 C eD-uise C 



A 

2^ 2^^<^ dzc 



S{ze,Ze)=0. (133) 



These are actually the cyclic orderings of the 3jm-symbols in the definition of the tetrahedral spin network function II59I I. 
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Proof of the proposition |4l We start from the proposition [3] and the equation (|116p . which we multiply by 
j^g ^'^"^."'"l'. to obtain the correct normalization in the sum over spins of the spinor network. We form the quantities 

1 le = l V^^e ) ■ 

^g^'^"!'"."'"!', 3'^^"''^''^'^ (ge) with the appropriate shifts on ji, j2, ie- All factors l/{2j)l can be absorbed this way. However, 
on the two terms which shift ji , Ilul'^t' + l)' leaves some factors of the type (j2 — 1/2+^6 — 1/2+^4 + 2). The latter can 
be generated from derivatives with respect to spinors, (2 + ^ X]eDD264 '^e^^'^)' since Eee — X^a ^e^z^ becomes in 
the spin representation E^e — 2je- The final term in (jll6p displays a dimension factor 2ji + 1, which can also be taken 
into account via the operator (1 + -^i^^) i^^ same with zi). Then, the proposition is obtained by summing 

over all spins (ji, . . . , jg). □ 
We have thus obtained the equation of the dynamics as a holomorphic differential equation. One can check explicitly 
that the Schwinger's generating function (I132p is a solution. It is however not clear to us whether other solutions can 
be found. In any case, looking for a solution as a power series will obviously reproduce in fine the recursion formula 
for the 6j-symbol. 

Notice that the generating function S is most naturally a function of the classical variables F*^, — [ze'\ze)- Hence, 
it is interesting to recast the above equation of the dynamics with derivatives with respect to F*^, . For instance, let 
us consider the explicit recursion given in appendix (|A10p . We multiply it with \/ (^2 + je + ji + 2) (j2 + je — ji + ^) 
and then by the appropriate factors to built from the three terms three generating functions. The spin dependent 
coefficients can be re-absorbed into multiplications and derivatives with respect to the F* variables. After some 
algebra, one gets 

Odd d d d 



dF^i dF,% dFi, ^''^'^ dF*^ dFi, OF, 



F21- 



65 



F*-r^ + F 



2^9^2*4 ^^9i^4*6 ^'9^6*2 



(^^a2|^ + l))5(^e.) = 0. (134) 



The first two terms of this equation can be understood as follows. Remember that a variable F*^, is associated to 
each couple of links which meet on a node of the graph. The closed expression for the generating function displays 
Ai, (I132p . the products of those variables which share a cycle, and a sum over all 3-cycles and 4-cycles. Besides, the 
operator which generates the equation we are looking at acts on a fixed cycle of the graph, chosen to be here (666162). 
The first term we observe (which comes from the shift of ji by +1/2) contains the derivatives with respect to the 
three variables which live on the chosen cycle, ^2*^, F*g, Fq2, which give a contribution 1 when acting on A4, 

d d d Odd 

dWidW^dFg^^dFi^dW^dFiy^' ^^^^ = ^ • ^^^^^ 

The second term is more interesting. Indeed, when the derivatives act on Ai, all cycles which do not contain the 
corresponding F*^, disappear. Hence, the derivatives with respect to ^"32 , -Fg 2 ) ^6*5 select one 4-cycle, that with the 
links (63626665), and then 

The idea is thus that this contribution transforms the 4-cycle into the 3-cycle (616563). 
We will discuss other equations satisfied by the generating function in the section IVBI 

V. COMPARISONS WITH OTHER EQUATIONS FOR THE BF PHASE 
A. Back to the plaquette operator, with spinors 



We have discussed in the section IIII Al the plaquette operator in terms of the character operator along a cycle, in 
the group and in the spin picture. Let us now explain how to get the analogous equation to (I55p in the basis of SU(2) 
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coherent spin networks s^^^''^"'^"^ (ge)- We work the character in the fundamental representation, from which all others 
can be obtained. It is then better to write the product of holonomies around p in terms of the bosonic operators, 

Xi(ffp) .^-'--^ = n s^'-'^^^y . (137) 

e— 1 -"^ 

As we have seen, the re-coupling process hi the spin basis enables to extract a 6j-symbol at each vertex. Hence, we 
will write the action of the character vertex by vertex, instead of a product over links. This gives 

Xi(5p)«^^-^-^=^= E n('^'^«^l^'='«^') .^s^^-^-^'K (138) 

Here we sum over all ee = 0, 1 to include all combinations of the operators E, F, at each vertex v. (e, e') denotes a 
pair of links on the boundary of p which meet at v, and so that e is outgoing and e' ingoing. In the usual spin basis 
we have already shown that the above operators produce the expected re-coupling coefhcients of the equation (I55p . 
To get the equation in the SU(2) coherent spin network basis from (1551) . one can use the following correspondence, 

(flelfle') 
[ae\ae') 

while ignoring the signs in (I55p . 
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(139) 
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je' + 2 je' ^ 2 
Je Je 2 ^1 



Je 



2 -^e 2 

Je - 2 ^1 



2je 



V2je' 2je {Ze'\Ze) , 



(141) 
(142) 



B. Different types of recursions and differential equations 

We have at our disposal the fundamental spin 1/2 recursion for the topological ground states, theorem [1] derived 
from a Hamiltonian -ff^i^p, and which we have shown to reproduce the standard flatness constraint in the section fill El 
This equation implies a difference equation on a single spin, which can be solved explicitly (cycle after cycle). So 
it amounts to completely solve the theory. Obviously the algebra generated by the new scalar Hamiltonian we have 
presented can be used to derive numerous recursion relations. 

Technically, the success of this Hamiltonian relies on the fact that the coefficients of the recursion, A±_o, are 
really non-trivial coefficients. Having such non-trivial coefficients is needed to relate three consecutive terms of the 
6j-symbols. But it could also be interesting to consider different recursions on the ground states which 

• involve numerous shifts on several spins and more terms than the usual three terms (and hence cannot be solved 
from a single initial condition), 

• but, in contrast with the usual recursions, have simple, or even trivial, coefficients, 

• and are generically much easier to derive than the fundamental spin 1/2 (or spin 1) recursion on the 6j-symbol. 

Such recursions should be considered analogous to Ward identities or Schwinger-Dyson equations, which encode very 
interesting properties of the theory and are simpler to deal with than the equations for the effective action itself. 

Since one can construct numerous recursion relations from successive applications of the Hamiltonian, we need some 
criteria to select the 'nice' and useful ones. We can distinguish the following cases 
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• the recursion is generated by a natural operator, like the plaquette operator (and from which a nice geometric 
interpretation can be found), 

• the recursion has trivial coefficients (this is only possible if the recursion contains a large number of terms) , 

• the recursion comes from a simple differential equation on the generating function. 

We have already encountered such recursions, that for the plaquette operators, interpreted as tent move evolutions 
in the section pil Ap . Clearly this operator is simpler that the new Hamiltonian, since this is just a SU(2) character, 
but it shifts all the spins around a face. We will now give other examples. Some of them were found in [H, [s^l- In 
particular, from an integral representation of the isosceles 6j-symbol, | a h | ; one can derive the following equation. 



Unlike the standard recurrence relations found in the literature, this one has trivial coefficients and all the spins are 
shifted. These features are achieved at the price of increasing the number of the terms in the relation, from the usual 
three or four to eight. To understand the status of this recursion, one should take the perspective of the spin 1/2 
recursions. Since the latter can be used to evaluate all 6j-symbols, the above equation should also be derived from 
them. That was done in [sO]- The main result of the derivation is that one has to i) take combinations of spin 
networks on which different Hamiltonians will act, ii) use spin 1/2 recursions on both ji, j2 successively, and iii) sum 
over different cycles and all values of the free parameters a, /3 in ((82|) . More precisely 



Q,0 = ±l/2 

is a recursion relation on the 6j-symbol, which shifts all the spins, and whose coefficients trivialize when choosing 
h = J5 = a, ja = je = so as to get (|143p. 

Another example was found in [48| from the integral representation of the squared 6j-symbol, and it was interpreted 
in the large spin limit as encoding the closure of the tetrahedron dual to the spin network graph. It also corresponds to 
a recursion on the amplitude of the model on the 3-sphere triangulated by two tetrahedra, in the spin representation 
(known as the spin foam representation of the partition function), and this is the first and only recursion of this type 
that we know so far. 

Finally, let us discuss the relationship between recursion relations in the spin network basis and differential equations 
on the associated generating function. The idea is simply that the coefficients of the recursion relation come from 
acting with derivatives on the generating function. The reason why the equation of the proposition |4] on the generating 
function is relatively complicated is that it encodes the fundamental recursion, so that several sums of partial derivatives 
are needed to create the non-trivial coefficients of the recursion. However, one can also observe that the generating 
function satisfies other equations which are simpler. 

As noticed in [t^, the special structure of the generating function, S — {\ + Al)^^ with M. being a sum over closed 
diagrams (here 3-cycles and 4-cycles) of the graph induces some interesting relations. Denote the operator 

Lee' = Ke , (145) 

and consider the fact that Lee' acts on by projecting it to only those cycles which contain F^^, , which we denote 
A'iee' = Lf.f.iM{F*). In the case of the 6j-symbol, each F*^, appears once in a 3-cycle and once in a 4-cycle (it can be 
read directly in (|132l) . Choose two couples of links (ab) and {ij) which meet on the same or on two different vertices. 
Then, 



A^ab Lij — Mij Lab 



S{F*) = . (146) 



This equation induces a recursion on the 6j-symbol, by expanding 5 as a power series and identifying the different 
orders. Its coefficients keep a simple form since we have only two first derivatives. However, the multiplication by Aiab 
and Mij produces numerous shifts on the 6j-symbols (and we prefer not to write it, the above form as a differential 
equation being much nicer). 



31 



C. Generalization of the new scalar Hamiltonians to arbitrary lattices with 3-valent vertices 

It is also worth comparing the case of the triangular lattice which we have now studied in details with the case of 
a generic lattice. The scalar Hamiltonians, with spin 1 or 1/2, we have introduced keep the same form when working 
with an arbitrary lattice. Denote the links on the boundary of of a plaquette p by (ei, . . . , e„) and consider any vertex 
on it, say v where ei, e„ meet, both ingoing, and take v as the reference vertex for the holonomy gp around p. Then, 
the spin 1 Hamiltonian is 

Hi]l = X„ • (id - Ad{gp)) = 1„ • Xi - X„ • Ad(5„_i ■ ■ ■ 92) X, , (147) 

where we have used the fact that X^, and Xg are related by parallel transport along e. Similarly, it is natural to extend 
the spin 1/2 Hamiltonian to 



^ y {z„\Zn){zi\zi} 

For a fixed p, important differences with the plaquette operator described in the previous section are 

• there is a single plaquette operator Xi(<?p)i but n Hamiltonians Hi]l, one for each node (this was a key point in 
the triangular case, to show that the Hamiltonians enable to get back to the flatness constraint in the section 
UlTEll . 

• each of them only affects (n — 2) spins, since the holonomies on the two edges meeting on the reference vertex 
do not enter the expression. 

Quite clearly, evaluating the actions of these operators on spin network states on the identity will produce recursion 
relations for Wigner 3nj-symbols, which are solutions to the topological phase. Recursion for the 9j and 12j-symbols 
are known from textbooks [5l|, and recursions for 15-symbols have been recently found in [sO]- They were discovered 
using an invariance under a four-dimensional Pachner move of 4d triangulations. We strongly think that the above 
Hamiltonians provide us with a good tool to generically understand the whole set of recursions for arbitrary Wigner 
symbols (this is work in progress (67|). Among the remaining issues is the fact that we only expect three constraints 
to be independent (since the curvature vanishing has three real components), and it is not clear why this is the case 
from the point of view of recursion relations. 

In the basis of coherent intertwiners used in the section IIVB| where there are no spins anymore, solutions to the 
topological phase should be given by generating functions of Wigner 3nj-symbols, i.e. in the Bargmann representation. 
These are holomorphic functions of 2L spinors, if L is the number of links. These generating functions are known from 
[73} (see also more recently 74]), but we have not analyzed the differential equations they satisfy so far. 



VI. CONCLUSION 



In this paper, we have given a very detailed account for the topological phase of the SU(2) BF model. While the 
continuum formulation is quite well- understood at a formal level (functional determinants in the path integral) , a full 
lattice description which can be explicitly handled was missing. It is important because the BF model is at the core 
of recent developments concerning topological aspects in condensed matter, quantum information and background 
independent quantization such as loop quantum gravity. 

Solutions are given by evaluations of spin networks, producing SU(2) Wigner coefficients. While that was expected, 
our approach, through the new Hamiltonian, goes deeper than previous works into the structure of the theory, by 
generating the basic recursions on Wigner 6j-symbols, completing this way the program of (3l| . 

In 2+1 dimensions, the holonomy- flux algebra representation of the new Hamiltonian emphasizes the geometric 
content of 2-1-1 gravity. As an outcome, we get a quantization of flat, Euclidean geometry in terms of recursions from 
group representation theory. We have given several kinds of such equations in the last section. Other equations which 
turn out to quantize simple geometric properties of flat space have appeared in [48l |50| . 
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As already noticed in [3l|, a surprising result is that there is not a single time reparametrization as expected 
from geometrodynamics, but instead different directions of evolutions which mix spatial diffeomorphisms and time 
reparametrization. 



In this topological model, the long-standing issue of the ambiguity on the spin of the Wilson loops used to regularize 
the curvature in the Hamiltonian gets a precise answer. The fundamental representation is the only one which enables 
to solve the model for the group SU(2) without having to supplement it with additional initial conditions by hand. 



The method we have used can be easily applied to other compact Lie groups. The cases of finite or Abelian Lie 
groups are even simpler. The point is that the Biedenharn-EUiott identity is a non-trivial relation which is part of 
the machinery of representation theory. It is always possible to get from it recursion relations like those we studied 
here, and they will always be generated by our new Hamiltonian. As we have argued, they are better suited than 
the recursions which come from the plaquette operators considered in [29], though part of the same mathematical 
framework. 



Our Hamiltonian can also be used on plaquettes with an arbitrary number of links, generating this way recursion 
relations which are satisfied by generic Wigner symbols. Some of them are known from some kind of Biedenharn-EUiott 
identity satisfied by the 9j-symbol [5l|, but the generic case is actually new. It would definitely be interesting to s tudy 
further those recursions, which certainly have a lot to do with integrability of Wigner coefficients as discussed in |43{ . 



When written with spinors, coherent states quantizations are available which change the recursions to partial 
derivative equations on holomorphic functions. Some of the most interesting features of those representations are the 
existence of a closed formula for the generating function of 6j-symbols (|132p and the simple equation (|146l) it satisfies. 
This is a special case of an analysis performed in [tsI for generic Wigner coefficients. The spinor variables have 
also been very useful in recent mathematics [s^, l42j |. We hope in the future to apply these tools in the context of 
topological order and the dynamics of loop quantum gravity. 



conditions. In the continuum, the BF model is known to feature exotic statistics 
dimension involve particles, strings and/or branes. The dynamics is presented ii 



75, 


76|,^ 


1 77 


, 78t 



which depending on the 
. It would certainly be 



interesting to introduce them in our framework, and to compare with lattice models exhibiting similar phenomena 
(topological order and exotic statistics involving branes) like [T^ where the model is exactly solvable. More generically, 
we think our model is robust enough to allow the study of theories expanded around the topological BF phase, as 
proposed in [s^l and realized in the continuum, Yang-Mills theory in [l9j . 
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Appendix A: Explicit formulae 

The explicit graspings 



^12 ijiisislO) = [ai\a2) tj.j^jM = \J {h + h + h + l) (ii + k - h) tfcifcsjslO) , (Al) 

Ei2inn]M = (ai|a2)tju2i3|0) = -5k^,n+k^k2,02~^ \J O'l " ^ + + l) ("Ji + ^ + Js) '^fcifcsialO) > (A2) 

^61 ij^JlJe = [oelSl) tj^jlje = --^fci Ji + i^e Je-i {h + k " je + l) {-jl + k + k) I'kijike , (A3) 

^16 h'Jiie = WISl] tj,j,j^ = S;^^j^_i5k^j^+i (- jl + k + k + 1) (il + k - k) l-kijike , (A4) 

^61 t-jViio = ("elSi) tj^j^j^ = S^^^ j^^iSk^j^^i (ii + k + k + 2) (jl - k + k + 1) ^fejiifce > (A5) 

EiQhtitJ<> = (ailae) tjj-j-je = (^fc^ i (5fe,,j,_ i ^ {h + k + Je + l) (ji - js + k) tfc^ijfee > (A6) 



P62hi3Sj4 = [S6|a2)tj*j|j4 = 5fej,+ i5fe6j6 + i (i2 + je + k + 2) (j2 + k - k + 1) ''k'^klu 

H&Hir.u = («6|S2]ij-j*,4 = -^;„j,_i<5;ej6-i ^ (i2 + k + k + 1) (i2 + k - k) i^iniu ■ 

The spin 1/2 recursions, 

^ + is + k) (ii + i5 - k + 1) (-J1 + i2 + k) {h - h + k + 1) j'^' ^ 



i2 - 5 J3 

is ie - 1 



+ V (ii + k - k) {-jl + k + k + 1) (ii + i2 - is) (ji + i2 + k + i) 



Jl - 5 i2 - ^ J3 

i4 is ie + 5 



(A7) 
(A8) 



rf..V02+-^6 + ^4 + i)(j2+j6-j4) I;;.; g| (A9) 
^ (ii - is + k) (ii + is + k + 1) (ii + i2 - is) (ii + i2 + k + 1) j"^' ^ ^ ^'^ ~ ^ ^.^•^^ |} " ° ' 



V (Jl - J5 + J6 + 1) (Jl + J5 + J6 + 2) (jl + 32 - JS + 1) (jl + J2 + JS + 2) j^^ + ^ + ^ ^.^^_[i| 

+ \/(i2 + k + k + 2) (i2 + k - k + 1) {^-^ ^-^l (AlO) 

- ^J{k + is - Je) + J5 + J6 + 1) (ji - 32 + Jz) (-ii + i2 + Jz + 1) l^'^r ^ ^2 + ^ ^^JS i| = . 

V (J'l - J5 + J6 + 1) (jl + J5 + J6 + 2) (-J1 + J2 + Js) {]l - 32 + JS + l) [^'^^ ^ ~ ^ ^.^-^^ij 

- '^ii \l{k-k + k){-k + k+k + i) {^-^ J (All) 
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Appendix B: Direct proof of the Wheeler-DeWitt equation 

We derive the difference equation (j82[) by direct evaluation in the spin network basis. Since i^J^g commutes with 
+ FeiEi2, we can use the resuh of the equation (|M|) . 
6 -, 

E tn djMUe)^ ^16 h2 Fl s\i;^ 

6 

be} e = l 

On the second step, we have just re-indexed the sum, by ji — > ji + 1/2, j2 — > j2 + 1/2,^6 je + 1/2, and we have 
used the following identities on the coefficients, 

At^^\ju.n - a, . . . , J6 - /?) = (-1)"+'^ J2, . . . , J6) , 

d,,+i + i J2 - a, J, -P) = d,, A^-r^-^\n,j2, . . . , Je) . 

Obviously, we have a similar equation for the action of -Fe i-E* 12-^26 ■ 

The key point is the action FjgF62 since we now have a different ordering. Basically, Fq2 will shift the spins j2,j6 

by +1/2 on the half- legs at the vertex v (figure [1]), and generate the coefficient ^q^'^'. Then, Fjg lowers these spins 
so that they get back to their initial value. The difference with the calculations of the proposition [5] is thus that the 
coefficient we will extract involves shifts of ^'2, je by -1-1/2 instead of —1/2. 

More precisely, we use the spin 1/2 graspings of the equations (|48| . (|47)) to get 

-F|ei?62 ^(51,.. .,56) = E^^-^^.2+irf.6rf,.+i [f[dj.] K^'+^O.)]' V'(je) 4;?t''"-''^^ • (B2) 

be} ^==1 

We can now sum the three contributions. They all have a common factor, dj^dj^j^i djgdjg_^_i [ne=i '^ie]^o^'^^(je), 
Hili'U ^(fi' ■ • ■ '56) = E dnd,,H ^^<^^seH \fl 4^'^^Oe) 4^1^"^ 

fie} e=l 

+ ,/. .1.1 ..1\ .f + . + l . . ^ , / . 1 



+ 2' J2 + 2' ■ • ■ ' J6 + 2) + ^L^i^'(il) ^(jl - 2' J2 + 2' • ■ • ' J6 + 2) + 4^'^^Jl) V-ljl, ■ • ■ , i6) 



The last step is to project the above equation on an arbitrary spin network state, using their natural inner product. 
It is well-known that they are orthogonal, and given their normalization (|60p. we get 

(4et^l i^i'lllL 1^) =4+'+^(jl) ^0-l,...,j6) 

+ ^{ji + 1,J2 + i, . . . , J6 + ^) + A^^^+\n) ^(ji -l,J, + ^,... ,J, + 1) . (B3) 

Appendix C: Additional representations 
1. Another generating function for 6j-symbols 



In this section, we use coherent intcrtwiners obtained from the coherent states for each Schwinger's boson. This 
enables to get rid of the spins, and to fully transform the recursion relation to a differential equation with respect to 
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spinors only. For a spinor z, we introduce on the space of irreducible SU(2) representations the following state 



E 



1 



E E 



|j, rn) . 



(CI) 



It is an eigenstate of a"^ with eigenvalue z^, and a"^^ acts like d^A as usual. Just like in the previsou section, coherent 
intertwiners are built by tensor products and group averaging of such states. Consider the node of the figure [T] where 
ei, 62, 63 meet and are all outgoing. The three half-lines carry spinors (zi, Z2, Z3), and the corresponding intertwiner is 



i{zi,Z2,zz)= ^ 



v/(2ji)!(2j2)!(2j3)! 



Pjlj2j3 (^1 J ^2 7 ^3) >'jlj2j3 



(C2) 



On a spin network graph, we associate one such intertwiner to each node, and glue them using the holonomies on 
the links. This is the natural gluing, where all irreducible representations are orthogonal. This leads to the notion of 
spinor networks, which we define for the moment as 



where the coherent spin network s^^ 



,ze,ze} jg defined in pi3p . Its evaluation on the unit gives 

6 



w = E [n 



1 



=1 (2je) 



3l J2 J3 

ii k k 



Proposition 5. The Hamiltonian equation H^^^'^^ ipS^^*"^^^ 



(5e 



^ '-A,B=±l/2 



dzf 



+ [Zl\ze) [Zl\z2) 



reads in this basis 





-)( 


-dz6 


dz2/ V 



s(Ze, Ze) = 



(C3) 



(C4) 



(C5) 



Proof of the proposition [5l This equation is easily obtained from the proposition |31 First correct the normaliza- 
tion to match that of the spinor symbol (IC4[) . This removes all spin dimension factors from the recursion (|116l) except 
(2ji -|- 1) on the term which has no shift. Nevertheless, this factor is simply generated by the operator (1 -|- En). In 



the coherent spinor basis, the latter is the operator {1 
the proposition. 



A -^1 Wz^^' then has just to sum over all spins to get 

□ 



2. U(N) coherent states 

The fact that the operators (-Eee') satisfy a u(iV) algebra naturally leads to the introduction of U(iV) coherent 
intertwiners [55| 

Jl+j2+j3 = Jl23eN yl\l^i{2je)l 

3 

= TH-^ I ^ (^(aelZe))''''" 0e |0)e ■ (C6) 

i2Ji23j! JsU(2) ^ 

Those vectors correspond to the decomposition of the space of intertwiners with N legs (above N = 3) Hn into sectors 
"Hj^"* in which the sum of all spins meeting at the node is fixed to je = J ^ N, 

H7v = 0Hj^^ with n^^^ = Inv[«)eHjJ . (C7) 
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Those U(A^) coherent states are eigenvectors of the total spin operator, 



(C8) 



and are covariant with respect to the natural action of SL(iV, (D) acting on the set of A'' spinors. 

Further, still assuming that all links are outgoing at v, we now that the operator i?ee' exchanges a quantum 1/2 
between the leg e and e', hence the total spin is kept fixed. The operator Fee' (respectively F^^,) destroys (respectively 
creates) a spin 1/2 on e and e', and thus changes the total spin by —1/2 (respectively +1/2), 



d 



Ei2 tji23(2;i, 2:2,23) ^^^2 tta '•^123(21,22,23) , 

A "^1 

-^^12 '-7123(21,22,23) = [21I22) tJi23-l(zi,Z2,Z3) , 



-^^12 '7123(21,22,23) = 



d 



dzi 



d 

dz2 



'./l23 + l(2l, 22, 23) 



(C9) 

(CIO) 
(Cll) 



One can obviously glue these intertwiners via holonomies, and get a U(A^) coherent 4J-symbol by setting the 
holonomies to the identity, 



^{,Jv^ 2e , Zg) 



1 



E 



1 



]1 J2 J3 



nnt 



\J„-2ie 



(C12) 



with the same notations as for the equation (|13ip . 

Proposition 6. The equation of the topological dynamics becomes in this basis 



d 



9zf 



'5( J123, "-^156, >^264 — 1, 'hib] 2e, ^e) 



■A,S=±l/2 

\zi\zq) [Zi\z2) 5(Ji23 - 1, Jl56 " 1,^264 " 1, J345 ; 2:e , 2e) 



r d 




-dze 


dz^ ( 



r(— ) 



The proof proceeds like in the previous sections, either by direct calculations from the Hamiltonian H^^^^^^^j^ 

terms of the operators E,F,F'', or starting from the same equation in another basis such as (|116l) . 

Let us describe the three terms of this equation. The first comes from that of the recursion formula for the 6j-symbol 
which shifts ji by +1/2, and j2,j6 by —1/2, so that at the level of the total spin on each vertex, only J264 is shifted. 
The second term is due to the term of the recursion which shifts ji, j2, all by ~l/2. This affects three total spins, 
•^123, </i56, t^264- Finally, the last contribution is that without any shifts on the 6j-symbol arguments. 

It is also possible to obtain the U(A^) coherent 4J-symbol from the generating function of 6j-symbols S{ze,Ze)- 
Indeed, introduce four real positive parameters (<^), one for each vertex, and substitute each spinor Ze with ^JT^Ze 
and Zg with ■y/tt(e)2e, where s(e), t(e) are the source and target vertices of the link e. Then, we get the desired function 
by differentiating iS(ze,Ze) (J^ + 1) times with respect to each before setting them to 0, 



S{Jy;Ze,Ze) =\Y[ 



(jy + ly.jj dti- 



S{\/ts{e) Ze, \/tI{I) Ze) 



(CM) 



To see that precisely, it is sufficient to focus on an intertwiner. From the definition (jC6|) and the fact that the SU(2) 
coherent intertwiner ij^{ze) is a homogeneous polynomial of degree 2(^g je) — 2Ji,, we get 



(,(v^zi, v^z2, v^zs) = "^{K + 1)! t'^ i_ft:(zi, Z2, Z3) , 



(C15) 
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so that 

1 d"^ 

(J+ 1)!J! dt' ^^2, Vizs) \^^^ = Z2, 23) • (C16) 

Finally notice that the generating functions we have observed so far are actually power series in the classical variables 
F*^, — [ze'\ze)- This strengthens the idea of js^ of quantizing the model with wave functions on the F* variables, i.e. 
representing the kinematical space of states on each iV-valent node "Hj^"* as 

nf> = {Pe F[F*], Vi e C P{tF*^,) = P{F*^,)} , (C17) 

the space of holomorphic, homogeneous polynomials of degree J in the F* variables. It could be interesting to study 
this proposal deeper, and to rewrite the Hamiltonian in this language, similarly to the equation (|134p obtained in 
the previous section on the generating function of 6j-symbols. However we will stop here the description of the 
Wheeler-DeWitt equation in different bases, and we will sketch different directions to generalize our present work. 



[1] X. G. Wen, "Topological orders and edge excitations in fractional quantum Hall states," Advances in Physics 44, 405 
(1995). 

[2] X. G. Wen and Q. Niu, "Ground-state degeneracy of the fractional quantum Hall states in the presence of a random 

potential and on high-genus Riemann surfaces," Phys. Rev. B 41, 9377 (1990). 
[3] J. Frohlich and T. Kerler, "Universality in quantum Hall systems," Nucl. Phys. B 354, 369 (1991). 

[4] X. G. Wen "Mean-field theory of spin-liquid states with finite energy gap and topological orders," Phys. Rev. B 44, 2664 
(1991). 

[5] V. Kalmeyer and R. B. Laughlin, "Equivalence of the resonating valence bond and fractional quantum Hall states," Phys. 
Rev. Lett. 59, 2095 (1987). 

[6] X. G. Wen, F. Wilczek and A. Zee, "Chiral spin states and superconductivity," Phys. Rev. B 39, 11413 (1989). 
[7] X. G. Wen, "Quantum orders and symmetric spin liquids," Phys. Rev. B 65, 165113 (2002). 
[8] M. Z. Hasan, and C. L. Kane, "Topological Insulators," Rev. Mod. Phys. 82, 3045 (2010). 

[9] G. Y. Cho and J. E. Moore, "Topological BF field theory description of topological insulators," Annals Phys. 326, 1515 
(2011) [arXiv:1011.3485 [cond-mat.str-el]]. 
[10] T. H. Hansson, V. Oganesyan and S. L. Sondhi, "Superconductors are topologically ordered," Annals of Physics 313, 497 
(2004). 

[11] M. C. Diamantini, P. Sodano and C. A. Trugenberger, "Superconductors with topological order," Eur. Phys. J. B 53, 19 

(2006) [arXiv:hep-th/0511192]. 
[12] G. T. Horowitz, "Exactly Soluble Diffeomorphism Invariant Theories," Commun. Math. Phys. 125, 417 (1989). 
[13] M. Blau and G. Thompson, "Topological Gauge Theories of Antisymmetric Tensor Fields," Annals Phys. 205, 130 (1991), 

M. Blau and G. Thompson, "A new class of topological field theories and the Ray-Singer torsion," Phys. Lett. B 228, 64 

(1989). 

[14] A. S. Cattaneo, P. Cotta-Ramusino, J. Frohlich and M. Martellini, "Topological BF theories in three-dimensions and 

four-dimensions," J. Math. Phys. 36, 6137 (1995) [arXiv:hep-th/9505027]. 
[15] N. Maggiore and S. P. Sorella, "Perturbation theory for antisymmetric tensor fields in four-dimensions," Int. J. Mod. Phys. 

A 8, 929 (1993) [arXiv:hep-th/9204044]. 
[16] F. Gieres, J. M. Grimstrup, H. Nieder, T. Pisar and M. Schweda, "Symmetries of topological field theories in the BV 

framework," Phys. Rev. D 66, 025027 (2002) [arXiv:hep-th/0111258]. 
[17] C. Lucchesi, O. Piguet and S. P. Sorella, "Renormalization and finiteness of topological BF theories," Nucl. Phys. B 395, 

325 (1993) [arXiv:hep-th/9208047]. 
[18] E. Witten, "On quantum gauge theories in two-dimensions," Commun. Math. Phys. 141, 153 (1991). 

[19] A. S. Cattaneo, P. Cotta-Ramusino, F. Fucito, M. Martellini, M. Rinaldi, A. Tanzini and M. Zeni, "Four-dimensional Yang- 
Mills theory as a deformation of topological BF theory," Commun. Math. Phys. 197, 571 (1998) [arXiv:hep-th/9705123]. 

[20] E. Witten, "(2-|-l)-Dimensional Gravity as an Exactly Soluble System," Nucl. Phys. B 311, 46 (1988). 

E. Witten, "Topology Changing Amplitudes in (2-|-l)-Dimensional Gravity," Nucl. Phys. B 323, 113 (1989). 

[21] R. De Pietri, L. Freidel, "so(4) Plebanski action and relativistic spin foam model," Class. Quant. Grav. 16, 2187-2196 
(1999). [gr-qc/9804071]. 

[22] L. Freidel, A. Starodubtsev, "Quantum gravity in terms of topological observables," [hep-th/0501191]. 



38 



[23] E. Dennis, A. Kitaev, A. Landahl and J. Preskill, "Topological quantum memory," J. Math. Phys. 43, 4452 (2002) 
[arXiv:quant-ph/0110143]. 

[24] A. Y. Kitaev, "Fault tolerant quantum computation by anyons," Annals Phys. 303, 2 (2003) [arXiv:quant-ph/9707021]. 

[25] J. C. Baez, "Spin network states in gauge theory," Adv. Math. 117, 253 (1996) [arXiv:gr-qc/9411007]. 

[26] J. C. Bacz, "An Introduction to spin foam models of quantum gravity and BF theory," Lect. Notes Phys. 543, 25 (2000) 

[arXiv:gr-qc/9905087]. 

[27] A. Perez, "Introduction to loop quantum gravity and spin foams," arXiv:gr-qc/0409061. 

[28] A. Ashtekar and J. Lewandowski, "Background independent quantum gravity: A Status report," Class. Quant. Grav. 21, 

R53 (2004) [arXiv:gr-qc/0404018]. 
[29] M. A. Levin and X. G. Wen, "String net condensation: A Physical mechanism for topological phases," Phys. Rev. B 71, 

045110 (2005) [arXiv:cond-mat/0404617]. 
[30] V. Aquilanti, A. C. P. Bitencourt, C. d. S. Ferreira, A. Marzuoli and M. Ragni, "Quantum and semiclassical spin networks: 

From atomic and molecular physics to quantum computing and gravity," Phys. Scripta 78, 058103 (2008) [arXiv:0901.1074 

[quant-ph]]. 

[31] V. Bonzom and L. Freidel, "The Hamiltonian constraint in 3d Riemannian loop quantum gravity," Class. Quant. Grav. 

28, 195006 (2011). [arXiv:1101.3524 [gr-qc]]. 
[32] L. Frcidcl and D. Louapre, "Ponzano-Roggo model revisited I: Gauge fixing, observables and interacting spinning particles," 

Class. Quant. Grav. 21, 5685 (2004) [arXiv:hep-th/0401076]. 
[33] J. W. Barrett and I. Naish-Guzman, "The Ponzano-Regge model," Class. Quant. Grav. 26, 155014 (2009) [arXiv:0803.3319 

[gr-qc]]- 

[34] B. Bahr, B. Dittrich and J. P. Ryan, "Spin foam models with finite groups," arXiv:1103.6264 [gr-qc]. 

[35] S. Garoufalidis, R. van der Veen and w. a. a. Zagier, "Asymptotics of classical spin networks," arXiv:0902.3113 [math.GT]. 
[36] F. Costantino and J. Marcho "Generating series and asymptotics of classical spin networks," arXiv:1103.5644 [math.GT]. 
[37] V. Aquilanti, H. M. Haggard, A. Hedeman, N. Jeevanjee, R. G. Littlejohn and L. Yu, "Semiclassical Mechanics of the 

Wigner 6j-Symbol," arXiv:1009.2811 [math-ph]. 
[38] J. Roberts, "Classical 6j-symbols and the tetrahedron," Geom. Topol. 3 (1999), 21-6 [arXiv:math-ph/9812013]. 
[39] V. Bonzom, E. R. Livine, M. Smerlak and S. Speziale, "Towards the graviton from spinfoams: The Complete perturbative 

expansion of the 3d toy model," Nucl. Phys. B 804, 507 (2008) [arXiv:0802.3983 [gr-qc]]. 
[40] M. Dupuis and E. R. Livine, "Pushing Further the Asymptotics of the 6j-symbol," Phys. Rev. D 80, 024035 (2009) 

[arXiv:0905.4188 [gr-qc]]. 

[41] J. W. Barrett, R. J. Dowdall, W. J. Fairbairn, H. Gomes, F. Hellmann and R. Pereira, "Asymptotics of 4d spin foam 

models," (2010). arXiv: 1003. 1886 [gr-qc]. 
[42] R. G. Littlejohn and L. Yu, "Semiclassical Analysis of the Wigner 9J-Symbol with Small and Large Angular Momenta," 

Phys. Rev. A 83, 052114 (2011) [arXiv: 1104. 1499 [math-ph]]. 

L. Yu, "Semiclassical Analysis of the Wigner 12J-Symbol with One Small Angular Momentum: Part I," arXiv:1104.3275 
[math-ph] . 

L. Yu, "Asymptotic Limits of the Wigner 15J-Symbol with Small Quantum Numbers," arXiv: 1104.3641 [math-ph]. 

[43] R. W. Anderson, V. Aquilanti and A. Marzuoli, "3nj Morphogenesis and Semiclassical Disentangling," J. Phys. Chem. A 
113 (2009) 15106. arXiv:1001.4386 [quant-ph]. 

[44] V. Bonzom, P. Fleury, "Asymptotics of Wigner 3nj-symbols with Small and Large Angular Momenta: An Elementary 
Method," [arXiv: 1108. 1569 [quant-ph]]. 

[45] R. J. Dowdall, H. Gomes and F. Hellmann, "Asymptotic analysis of the Ponzano-Regge model for handlebodies," J. Phys. 
A 43, 115203 (2010) [arXiv:0909.2027 [gr-qc]]. 

[46] K. Schulten and R. G. Gordon, "SEMICLASSICAL APPROXIMATIONS TO 3J AND 6J COEFFICIENTS FOR QUAN- 
TUM MECHANICAL COUPLING OF ANGULAR MOMENTA," J. Math. Phys. 16, 1971 (1975). 

[47] M. Dupuis and E. R. Livine, "The 6j-symbol: Recursion, Correlations and Asymptotics," Class. Quant. Grav. 27, 135003 

(2010) [arXiv:0910.2425 [gr-qc]]. 

[48] V. Bonzom and E. R. Livine, "Yet Another Recursion Relation for the 6j-Symbol," arXiv:1103.3415 [gr-qc]. 

[49] V. Bonzom, "Spin foam models and the Wheeler-DeWitt equation for the quantum 4-simplex," Phys. Rev. D84, 024009 

(2011) . [arXiv:1101.1615 [gr-qc]]. 

[50] V. Bonzom, E. R. Livine and S. Speziale, "Recurrence relations for spin foam vertices," Class. Quant. Grav. 27, 125002 

(2010) [arXiv:0911.2204 [gr-qc]]. 
[51] D. A. Varshalovich, A. N. Moskalcv and V. K. Khersonsky, "QUANTUM THEORY OF ANGULAR MOMENTUM: 

IRREDUCIBLE TENSORS, SPHERICAL HARMONICS, VECTOR COUPLING COEFFICIENTS, 3NJ SYMBOLS," 

SINGAPORE, SINGAPORE: WORLD SCIENTIFIC (1988) Slip 
[52] K. Noui and A. Perez, "Three-dimensional loop quantum gravity: Physical scalar product and spin foam models," Class. 

Quant. Grav. 22, 1739 (2005) [arXiv:gr-qc/0402110]. 
[53] L. Smolin, "The Classical limit and the form of the Hamiltonian constraint in nonperturbative quantum general relativity," 



39 



arXiv:gr-qc/9609034. 

[54] L. Preidel and E. R. Livine, "The Fine Structure of SU(2) Intertwiners from U(N) Representations," J. Math. Phys. 51, 

082502 (2010) [arXiv:0911.3553 [gr-qc]]. 
[55] L. Freidel and E. R. Livine, "U(N) Coherent States for Loop Quantum Gravity," J. Math. Phys. 52, 052502 (2011) 

[arXiv: 1005.2090 [gr-qc]]. 

[56] E. F. Borja, L. Freidel, L Garay and E. R. Livine, "U(N) tools for Loop Quantum Gravity: The Return of the Spinor," 

Class. Quant. Grav. 28, 055005 (2011) [arXiv: 1010.5451 [gr-qc]]. 
[57] L. Freidel and S. Speziale, "From twistors to twisted geometries," Phys. Rev. D 82, 084041 (2010) [arXiv:1006.0199 [gr-qc]]. 
[58] E. R. Livine, J. Tambornino, "Loop gravity in terms of spinors," [arXiv: 1109.3572 [gr-qc]]. 

[59] E. R. Livine, S. Speziale, J. Tambornino, "Twistor Networks and Covariant Twisted Geometries," [arXiv:1108.0369 [gr-qc]]. 

[60] E. R. Livine, J. Tambornino, "Spinor Representation for Loop Quantum Gravity," [arXiv:1105.3385 [gr-qc]]. 

[61] H. Ooguri, "Partition functions and topology changing amplitudes in the 3-D lattice gravity of Ponzano and Regge," Nucl. 

Phys. B 382, 276 (1992) [arXiv:hep-th/91 12072]. 
[62] V. Bonzom and M. Smerlak, "Bubble divergences from twisted cohomology," arXiv:1008.1476 [math-phj. 
[63] L. Freidel and S. Speziale, "Twisted geometries: A geometric parametrisation of SU(2) phase space," Phys. Rev. D 82, 

084040 (2010) [arXiv:1001.2748 [gr-qc]]. 
[64] V. Bonzom, "Geometric quantique dans les mousses de spins: De la theorie topologique BF vers la relativite generale," 

arXiv: 1009.5100 [gr-qc]. 

[65] B. Bahr and B. Dittrich, "(Broken) Gauge Symmetries and Constraints in Regge Calculus," Class. Quant. Grav. 26, 225011 

(2009) [arXiv:0905.1670 [gr-qc]]. 
[66] B. Dittrich and J. P. Ryan, "Phase space descriptions for simplicial 4d geometries," Class. Quant. Grav. 28, 065006 (2011) 

[arXiv:0807.2806 [gr-qc]]. 

[67] V. Bonzom and P. Floury, "A geometric approach to the evaluation of classical spin networks," in preparation. 
[68] E. R. Livine and S. Speziale, "A New spinfoam vertex for quantum gravity," Phys. Rev. D 76, 084028 (2007) 
[arXiv:0705.0674 [gr-qc]]. 

[69] F. Conrady and L. Freidel, "Quantum geometry from phase space reduction," J. Math. Phys. 50, 123510 (2009) 
[arXiv:0902.0351 [gr-qc]]. 

[70] L. Freidel, K. Krasnov and E. R. Livine, "Holomorphic Factorization for a Quantum Tetrahedron," Commun. Math. Phys. 

297, 45 (2010) [arXiv:0905.3627 [hcp-th]]. 
[71] A.C.T. Wu, "Structure of the Wigncr 9j Coefficients in the Bargmann Approach," J. Math. Phys. 13, 84-90 (1972). 
[72] J. Schwingcr, "On angular momentum," Report US AEC NYO-3071 (1952), in Quantum Theory of Angular Momentum, 

eds. LC Biedenharn and H. van Dam (New York: Academic, 1965). 
[73] J. J. Labarthe, "Generating Functions for the Coupling Recoupling Coefficients of SU(2)," J. Phys. A 8, 1543 (1975). 
[74] O. Schnctz, "Generating Functions for Multi-j-Symbols," [arXiv:math-ph/9805027]. 

[75] M. Bergeron, G. W. Scmenoff and R. J. Szabo, "Canonical bf type topological field theory and fractional statistics of 

strings," Nucl. Phys. B 437, 695 (1995) [arXiv:hep-th/9407020]. 
[76] J. C. Baoz, D. K. Wise and A. S. Crans, "Exotic statistics for strings in 4d BF theory," Adv. Theor. Math. Phys. 11, 707 

(2007) [arXiv:gr-qc/0603085]. 

[77] J. C. Bacz and A. Perez, "Quantization of strings and branes coupled to BF theory," Adv. Theor. Math. Phys. 11, 3 (2007) 
[arXiv:gr-qc/0605087]. 

[78] W. J. Fairbairn, A. Perez, "Extended matter coupled to BF theory," Phys. Rev. D78, 024013 (2008). [arXiv:0709.4235 

[gr-qc]]. 

[79] H. Bombin and M. A. Martin-Delgado, "Exact Topological Quantum Order in D=3 and Beyond: Branyons and Brane-Net 

Condensates," Phys. Rev. B 75, 075103 (2007) [arXiv:cond-mat/0607736]. 
[80] L. Freidel and K. Krasnov, "Spin foam models and the classical action principle," Adv. Theor. Math. Phys. 2, 1183 (1999) 

[arXiv:hep-th/9807092] . 



